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1.  SUMMARY  OF  WORK 


Many  of  the  past  year's  research  results  have  been  described 
in  seven  publications.  Reprints  of  those  publications  are  included  as 
Appendix  A of  this  report,  and  will  be  referenced  as  A.1-A.7  in  the 
following. 

The  specific  topics  addressed  during  the  past  year  are: 

1.  Effects  of  alloy  scattering  on  transport  and  device  characteristics 
[A.l]. 

2.  Origins  of  negative  resistance  in  polar  III-V  semiconductors  [A. 2] 

*3.  Transport  in  Ga^  In  P,  As  [A. 3]. 

^ 1-x  X 1-y  y 

*4.  Qualitative  analysis  of  thin  GaN  films  [A. 4] 

5.  Velocity-field  characteristics  of  GaAs  with  the  F-L-X  conduction 
band  ordering  recently  demonstrated  by  Aspnes"*"  [A.  5] 

6.  Energy  bandgap  and  lattice  constant  contours  of  III-V  quaternary 

alloys.  The  contours  for  the  2:2  alloy  systems  of  the  form 

A,  B C,  D were  calculated  first,  and  were  published  this  year 
1-x  X 1-y  y 

[A. 6].  The  contours  for  the  1:3  and  3:1  alloys  were  calculated 
and  have  been  sxibmitted  for  publication. 

7.  Performance  of  microwave  MESFET's  constructed  from  III-V  compounds 
and  alloys  [A. 7]. 

8.  Refinement  and  extension  of  the  Monte-Carlo  transport  programs 

(a)  Inclusion  of  p-state  mixing 

(b)  Work  toward  modeling  of  degenerate  materials 

(c)  Transient  field  program 

(d)  Position-dependent  field  program 

9.  Systematic  search  for  promising  materials  (calculation  of  v-E 
characteristics  of  III-V  compounds  and  alloys). 

*This  work  was  done  during  the  grant  period  1 Jan  76-31-Dec  76. The 
papers  were  published  this  year  and  are  included  for  completeness. 

+D.  E.  Aspnes,  Phys.  Rev.  B 5331  (1976). 


The  past  year's  research  efforts  have  been  quite  fruitful.  Considerable 
progress  has  been  made  in  the  areas  of  high-field  transport  and  the 
identification  of  materials  that  show  promise  for  microwave  device  application. 
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Apj)endix  A:  Publications  for  1977 


The  following  publications  are  collected  in  this  Appendix  for  ready  reference. 

1.  Littlejohn,  M.  A.,  J.  R.  Hauser,  T.  H.  Glisson,  D.  K.  Ferry,  and 
J.  W.  Harrison.  "Alloy  Scattering  and  High-Field  Transport  in 
Ternary  and  Quaternary  III-V  semiconductors".  Paper  presented  at 
the  1st  Int.  Conf.  on  Hot  Electrons,  Denton,  EX,  July  1977. 

Also  accepted  for  publication  in  Dec.  77  issue  of  Sol.  St.  Electronics. 

2.  Hauser,  J.  R.  , M.  A.  Littlejohn,  and  T.  H.  Glisson.  "Negative 
Resistance  and  Peak  Velocity  in  the  Central  (000)  Valley  of  III-V 
Semiconductors".  Paper  presented  at  the  1st  Int.  Conf.  on  Hot 
Electrons,  Denton,  TX,  July  1977. 


Littlejohn,  M.  A.,  J.  R.  Hauser,  and  T.  H.  Glisson.  "Velocity-Field 
Characteristics  of  Ga,  In  P,  As  Quaternary  Alloys".  Appl.  Phys. 
L.,  v30,  no.  5 (March^l^VV^ . ^ 


4.  Andrews,  J.  E.,  A.  P.  Duhamel,  and  M.  A.  Littlejohn.  "Qualitative 
Analysis  of  Thin  Gallium  Nitride  Films  with  Secondary  Ion  Mass 
Spectrometry".  Anal.  Chem.,  v49 , no.  11  (September  1977). 

5.  Littlejohn,  M.  A.,  J.  R.  Hauser,  and  T.  H.  Glisson.  "Velocity- 
Field  Characteristics  of  GaAs  with  r^L?-X^  Conduction-Band  Ordering". 
J.  Appl.  Phys.,  v48,  no.  II  (November  toVv). 

6.  Glisson,  T.  H.,  J.  K.  Hauser,  and  M.  A.  Littlejohn.  "Energy  Bandgap 
and  Lattice  Constant  Contours  of  III-V  Quaternary  Alloys".  Accepted 
for  publication  in  Dec.  1977  issue  of  J.  Elec.  Mat. 

7.  Littlejohn,  M.  A.,  J.  R.  Hauser,  T.  H.  Glisson,  and  L.  A.  Arledge. 
"High  Field  Transport  and  Microwave  MESFET  Performance  in  Ternary 
and  Quaternary  III-V  Semiconductors",  Accepted  for  publication  in 
the  Proceedings  of  the  Sixth  Biennial  Cornell  Microwave  Conference. 
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ALLOY  SCATTERING  AND  HIGH  FIELD  TRANSPORT 
IN  TERNARY  AND  QUATERNARY  III-V  SEMICONDUCTORS* 

M.  A.  Littlejohn,  J.  R.  Hauser,  T.  H.  Glisson 
Electrical  Engineering  Department 
N.  C.  State  University 
Raleigh,  NC  27607 

D . K . F erry 

Office  of  Naval  Research 
Arlington,  VA  22217 

J.  W.  Harrison 
Research  Triangle  Institute 
Resecurch  Triangle  Park,  NC  27709 


ABSTRACT 


A technique  is  described  for  the  estimation  of  the  influence  of 
random  potential  alloy  scattering  on  the  high  field  transport  properties 
of  quaternary  III-V  semiconductors  obtained  by  Monte  Carlo  simulation. 

The  approach  is  based  on  an  extension  of  a theoretical  model  for  scattering 
in  the  ternary  alloys.  The  magnitude  of  the  scattering  potential  is  an 
important  parameter  in  alloy  scattering,  and  three  proposed  models  for 
calculating  this  potential  are  discussed.  These  are  the  energy  bandgap 
difference,  the  electron  affinity  difference,  and  the  heteropolar  energy 
difference  for  the  appropriate  binary  compounds . 

The  technique  is  used  in  the  Monte  Carlo  method  to  study  the  influence 
of  alloy  scattering  on  the  transport  properties  of  III-V  quaternary  alloys. 
The  results  of  this  study  are  used  in  a device  model  to  estimate  device 
parameters  for  FET's. 


*ThiE  work  was  supported  by  research  contract  No.  N00014-70-A-0120-004  from 
the  Office  of  Naval  Research,  Arlington,  VA,  and  by  research  contract 
No.  F3361h-76-C-1265  from  the  Air  Force  Avionics  Laboratory,  Wright  Patterson 
AfB,  Ohio. 


INTRODUCTION 


In  a semiconductor  solid  solution  alloy,  the  scattering  of  free 
carriers  due  to  the  deviations  from  the  perfect  periodicity  of  the 
virtual  crystal  model,  originally  conceptualized  by  Nordheim  [1],  can 
be  called  random  potential  alloy  scattering.  The  discussion  of  alloy 
scattering  in  relation  to  experimental  electron  mobility  in  III-V  semi- 
conductor ternary  alloys  has  been  based  on  an  unpublished  result  of 
Brooks  [2],  as  quoted  and  used  for  example  in  the  results  of  Tietjen  and 
Weisberg  [3],  Makowski  cUid  Glicksman  [43,  Glicksman,  et  al.  [5],  and 


Takeda,  et  al.  [6].  Recently,  theoretical  calculations  [7,8]  have 

elaborated  on  physical  models  for  alloy  scattering  in  ternary  III-V 

semiconductors,  and  these  models  have  been  applied  in  the  Monte  Carlo 

method  to  include  the  effects  of  alloy  scattering  in  high  field  transport 

calculations  for  InP,  As  ternary  alloys  [9]  and  Ga  In  As 

X .L“X  X X y y 

quaternary  alloys  [10]. 

The  III-V  quaternary  semiconductor  alloys  are  becoming  of  great 
technological  importance,  and  offer  for  device  applications  the  unique 
feature  that  the  energy  band  gap  can  be  varied  while  maintaining  a fixed 
lattice  constant  [11]  by  varying  the  alloy  composition.  It  is  important 
in  these  materials  to  be  able  to  estimate  the  influence  of  alloy  scattering 
on  the  electronic  transport  properties.  The  Monte  Carlo  method  [12]  is 
one  of  the  more  reliable  techniques  for  nviking  such  an  estimate,  and  this 
technique  does  become  a predictive  tool  for  III-V  ternary  and  quaternary 
materials  as  long  as  reliable  material  property  data  are  available  for  the 
binary  constituents  [7-lOj. 


The  purpose  of  this  paper  is  to  present  a technique  for  estimating 
the  effect  of  alloy  scattering  on  transport  properties  of  III-V  ternary 
and  quaternary  semiconductors.  This  procedure  has  been  used  previously 
[9,10],  and  will  be  discussed  in  detail  in  this  paper.  The  results 
obtained  from  the  method  applied  in  Monte  Carlo  transport  calculations 
are  also  used  to  predict  the  influence  of  alloy  scattering  on  device 
performance.  The  intent  is  to  establish  some  reasonable  upper-  and 
lower-bounds  for  the  effect  of  this  scattering  process  in  the  ternary 
and  quaternary  materials. 
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TERNARY  ALLOY  SCATTERING  POTENTIAL 

The  electron  scattering  rate  in  ternary  semiconductors  due  to 
alloy  scattering  is  important  for  the  development  of  the  quaternary 
scattering  rate.  Of  particular  importance  is  the  scattering  potential 
used  in  this  calculation.  The  scattering  rate  (transition  probability) 
for  a ternary  alloy  (A^^  with  nonparabolic  energy  bands  is  given  in 

r.m.k.s.  units  by  [7-9]. 

3/2 

— Ud-x):  Y(e)  a|aul^S(a)  1) 

"ta  e/5  dc  ' ' 

where  the  notation  used  is  described  in  Reference  10.  Here  S(a)  is  an 

energy-dependent  parameter  which  describes  the  effect  of  ordering  on  the 

scattering  rate  [13].  In  general,  0 < S <1,  where  S=0  for  a perfectly 

ordered  ternary  (superlattice)  and  S=1  for  a completely  random  alloy. 

The  evaluation  of  S is  complex  and  depends  on  the  scattering  potential, 

AU.  In  this  paper  the  two  cases  S=0,  1,  along  with  different  scattering 

potentials  for  the  case  S=l,  will  be  used  to  give  lower-  euid  upper-bound 

effects  of  alloy  scattering,  according  to  results  obtained  from  the  Monte 

Carlo  method. 

Besides  the  effect  of  ordering,  the  most  significant  parameter  in 
Eqn.  1 is  the  scattering  potential,  AU.  The  derivation  of  Eqn.  1 
is  based  on  the  Mott-inner  potential  [8],  since  this  potential  results 
in  a relaxation  time  which  leads  to  the  accepted  temperature  dependence  of 
electron  mobility  due  to  alloy  scattering  [5,7].  Previous  calculations  have 
used  either  a)  the  difference  in  energy  band  gaps  between  the  binary 


constituents  [2-6],  or  b)  the  difference  in  electron  affinities  between  the 

binary  constitutents  [7-10]  for  this  scattering  potential. 

Recently,  Ferry  [14]  has  suggested  another  form  for  the  scattering 

potential  based  on  the  electronegativity  theory  of  Phillips  [15],  In  this 

work,  the  energy  difference  between  the  bonding  and  antibonding  hybridized 
3 

molecular  (sp  ) orbitals  of  a tetrahedrally  coordinated  crystal  can  be 
decomposed  into  contributions  due  to  symmetric  and  antisymmetric  potentials 
within  a unit  cell.  This  relation  is  described  by  [15] 


2 2 2 


(2) 


where  E is  the  energy  difference  between  the  bonding  and  antibonding 
G 


molecular  states,  E is  the  homopolar  (symmetric)  part  and  C is  the  heteropolar 

H 


(antisymmetric)  part  of  this  energy  difference.  The  heteropolar  energy  C 
represents  the  charge  transfer  or  ionic  contribution  to  E^,  and  if  in  the 
alloy  all  bond  lengths  are  equal,  then  it  can  be  assumed  that  any  fluctuations 
in  the  crystal  potential  arise  from  fluctuations  in  C [16],  This  approach 
has  had  some  limited  success  in  explaining  energy  bandgap  "bowing"  effects 
in  ternary  semiconductors  [16,17]. 

In  Appendix  A an  expression  for  the  scattering  potential  in  a ternary 


alloy  based  on  the  electronegativity  theory  is  given  [14].  This 


expression  is 


AU 


= [ — ]exp[-k  R ] 

EN  47rE  r,  r„  ^ s A 

o A B 


(3) 


Each  quantity  in  this  expression  is  defined  in  Appendix  A. 


Table  1 lists  values  of  AU„.,  for  several  ternary  III-V  alloys,  along 

EN 


with  the  energy  gap  difference,  and  AU^^  , the  electron  affinity 


EA’ 


difference.  The  values  of  AU^^  are  given  in  Table  1 for  x=0.5,  since 


there  is  a slight  functional  dependence  of  on  the  alloy  composition  x. 

This  is  illustrated  in  Figure  1 where  the  compositional  dependence  of 


is  shown  for  several  ternary  alloys. 
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Figure  2 shows  the  effect  of  the  use  of  these  three  scattering 

potentials  on  the  velocity  field  characteristics  for  Ga^  ^As  obtained 

from  Monte  Carlo  calculations.  Also  shown  in  this  figure  is  the 

characteristic  for  S=0  (no  alloy  scattering),  and  the  curve  for  GaAs 

for  comparison.  These  calculations  have  used  the  F-L-X  conduction  band 

ordering  proposed  for  GaAs  by  Aspnes  [18],  and  GaAs  material  parameters 

resulting  in  a good  fit  to  experimental  data  [19].  Also  given  in 

Figure  2 is  a tabulation  of  the  low-field  electron  mobility  for  each 

scattering  potential.  In  the  recent  work  of  Takeda,  et  al.  [6]  the 

2 

experimental  Hall  mobility  for  Ga  ^^In  is  8450  cm  /volt  sec. 

The  Monte  Carlo  drift  mobilities  calculated  by  the  maximum  likelihood 

2 2 

estimation  technique  [9,10]  are  7350  cm  /volt  sec,  8300  cm  /volt  sec,  and 
2 

8900  cm  /volt  sec  using  I'espectively , in  the 

calculations.  In  general,  the  Hall  mobility  is  greater  than  the  drift  mobility, 
so  that  based  on  these  low  field  mobility  calculations  the  scattering  potential 
is  possibly  closer  to  that  predicted  by  the  electron  affinity  difference  for 
this  material. 

If  one  examines  the  trends  for  the  eighteen  possible  III-V  ternary 
alloys  it  appears  that  often  < AU^.^.  However,  this  is  not 

always  true,  as  can  be  seen  by  the  examples  chosen  for  Table  I.  At  best, 
this  ordering  of  the  AU's  seems  fortuitous,  and  the  ordering  between 

and  AU^,^  is  especially  in  question.  The  band  gap  differences  are  the 

i 

most  accurately  known  parameters,  while  there  are  uncertainties  in  the  | 

I 

available  values  of  electron  affinities  [20]  and  covalent  radii  [21].  Thus  I 


the  ordering  in  Table  1 could  be  the  result  of  experimental  variations 
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especially  between  the  and  the  At  the  present  time  there  are 

no  firm  theoretical  or  experimental  reasons  for  choosing  either  scattering 
potential  to  evaluate  the  scattering  rate  due  to  alloy  scattering  in  the 
ternaries.  The  point  to  be  made  is  that  alloy  scattering  as  used  in  the 
Monte  CcUT'lo  method  and  based  on  Eqn.  1 has  a very  detrimental  effect  on 
the  transport  properties  of  III-V  ternary  alloys  and  thus  will  be  a factor 
in  their  use  in  devices,  if  the  proposed  models  are  correct. 


QUATERNARY  ALLOY  SCATTERING  RATE 


Appendix  B develops  an  extension  of  the  ternary  alloy  scattering 
rate  model  to  a quaternary  material,  ^^-x^x^l-y^y ’ ^ ® 

group  III  atoms  and  C and  D are  group  V atoms.  This  relation,  which 
applies  to  the  case  where  the  A and  B atoms  are  randomly  distributed 
on  the  group  III  sites  and  the  C and  D atoms  are  randomly  distributed 
on  the  group  V sites,  is  given  by 


where 


^ = k|aU  (x,y)|^, 

QA  ^ 


sn  -fi 


(4) 


and 


|A'J^(x,y)|^  = x(l-x)y^|AU^gQ|^  + x(l-x)(l-y)^  |AU^g^  | 


+ y(l-y)x^|AUgj,p|  + y(l-y)(l-x)^|AU^ 


ACD' 


Here  the  effective  mass,  m*,  and  the  primitive  cell  volume,  C,  are  calculated 

according  to  an  interpolation  procedure  described  previously  [10], 

Table  2 lists  the  quaternary  alloy  scattering  potential  AU^^  of  Eqn.  4 

at  the  mid- composition  range  x=y=0.5  for  three  quaternary  alloys,  Ga^_^In^P^_^ 

Ga,  In  P,  Sb  , and  Al,  In  P,  As  for  each  of  the  three  ternary  scattering 
1-x  X 1-y  y 1-x  X 1-y  y ^ ^ 

potentials  described  in  the  last  section.  When  comparing  the  results  of 
Tables  1 and  2,  the  numbers  given  in  Table  1 should  be  divided  by  1/4  since 
the  factor  x(l-x)  in  Eqn.  1 is  not  included  in  the  Table  1 entries,  whereas 
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the  analogous  factor  is  included  in  the  entries  of  Tiible  2.  In  general, 
the  quaternary  scattering  potential  (not  the  scattering  rate)  is  larger 
than  the  scattering  potential  for  any  of  the  four  ternaries  of  which  the 
quaternary  is  composed. 

Figure  3 shows  a plot  of  the  quaternary  alloy  scattering  parameter  surface 

2 

(AU„,)  for  Ga,  In  P,  As  in  the  compositional  plane  (0  < x < 1,  0 < y < 1) 

QA  1-x  X 1-y  y _ _ _ _ 

for  the  case  where  the  ternary  scattering  potentials  are  taken  as  the 

electron  affinity  differences.  This  figure  shows  a relative  minimum  of  the 

2 

quaternary  alloy  scattering  parameter,  AU^. , along  a region  which  is  very 

yn 

close  to  the  compositions  required  for  lattice  matching  this  quaternary 
to  InP  substrates.  This  contributes  to  the  large  predicted  peak  velocity 
[10 J of  this  particular  quaternary  alloy  matched  to  InP.  However,  if  the 
energy  gap  difference  is  used  for  the  scattering  potential  the  shape  of  the 
curve  is  shifted  to  a less  favorable  situation  for  minimum  alloy  scattering 
using  In?  as  a substrate.  Also,  for  cornpositions  away  from  this  region 
of  minimum  alloy  scattering,  the  effects  of  alloy  scattering  are  more 
detrimental.  This  can  be  seen  in  Figure  4,  where  the  velocity-field 
characteristic  for  Ga  ccin  P As  are  shown.  Without  alloy  scattering 

• DO  aHO  aJ. 

these  characteristics  are  nearly  as  attractive  for  device  applications  as 
those  previously  reported  for  Ga  ^.yln  ^As  However,  alloy  scattering 
has  a much  more  detrimental  effect  on  the  velocity-field  curves  for  this 
particular  composition.  Also,  for  the  composition  shown  in  Figure  4,  the 
effect  of  alloy  scattering  due  to  the  use  of  AU^^^  and  AU^,^  are  almost 
identical,  and  only  one  velocity-field  curve  is  shown  for  the  calculations 
made  from  each  of  these  alloy  scxittcring  potentials. 
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Other  different  effects  of  alloy  scattering  in  quaternary  alloys 
are  illustrated  in  the  velocity-field  characteristics  shown  in  Figures 

5 and  6 for  A1  25In.75P.25''^  75  25^".  75^  84^^16’  respectively. 

Since  our  interest  in  these  materials  has  been  primarily  for  FET's, 

17  -3 

the  velocity-field  curves  are  given  for  a doping  level  of  10  cm 
It  is  seen  that  when  alloy  scattering  is  not  used  in  the  Monte  Carlo 
simulations,  the  general  features  of  these  velocity-field  curves  eire  most 
attractive  for  a wide  veiriety  of  solid  state  Gunn-effect  electronic 
devices  [24],  These  features  include  large  low-field  mobility,  high  peak 
velocity,  low  threshold  field,  large  energy  band  gap,  large  intervalley 
energy  band  separation,  large  negative  differential  mobility,  and  large 
peak-to- valley  drift  velocity  ratio.  When  alloy  scattering  is  included, 
its  effects  are  generally  detrimental  to  all  these  desirable  properties. 

The  extent  of  this  degradation  is  certainly  open  to  question,  although 
the  use  of  Monte  Carlo  method  with  the  present  alloy  scattering  model  can 
provide  a very  good  estimation  of  the  range  of  this  degradation. 

At  the  same  time  the  results  of  the  Monte  Carlo  analysis  can  be  used 
in  device  models  to  estimate  device  f igures-of-merit  and  to  examine  their 
degradation  from  alloy  scattering.  This  is  illustrated  in  Table  3,  where 
the  FET  model  of  Lehovec  and  Zulceg  [25]  has  been  used  to  calculate  some 
device  parameters  for  MESFET's  using  the  ternary  and  quaternary  materials 
discussed  in  this  paper.  While  this  particular  device  model  includes 
velocity  saturation,  but  does  not  include  a negative  differential  mobility, 
it  has  recently  been  used  as  a design  model  for  GaAs  MESFET’s  [26],  with 
good  results.  The  important  comparison  to  be  made  is  that  between  the 
device  parameters  for  GaAs  and  the  ternary  and  quaternary  materials,  and 
to  consider  the  effect  of  alloy  scattering  on  these  device  figures-of-merit . 
These  three  materials  have  properties  which  suggest  possible  improvements 


I 


in  MESFET  performance  over  C.aAs  (although  the  bandgap  of  0.8  eV  for 


Ga  cln  may  be  somewhat  too  low),  and  these  first-order  model  calculations 

• D • D 

verify  this  suggestion.  Depending  on  the  amount  of  alloy  scattering  and  the 


physical  correctness  of  the  proposed  model,  alloy  scattering  reduces  the 
advantages  the  materials  offer.  However,  even  considering  the  uncertainty 

.1 

in  the  magnitude  of  the  alloy  scattering  the  ternary  and  quaternary  device 
parameters,  such  as  f_  in  Table  3,  are  seen  to  be  improved  considerably 

i I 

I 

over  the  GaAs  device  parameters. 
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This  paper  has  presented  techniques  which  can  be  used  in  the  Monte  Carlo  | 

method  to  estimate  the  effects  of  random  potential  alloy  scattering  on  the  : 

high  field  transport  properties  of  ternary  and  quaternary  III-V  semiconductors. 

These  materials  have  many  properties  which  suggest  that  their  utilization  in 
Gunn-effect  electron  devices  can  improve  presently  achievable  device 
performance.  At  the  same  tir-',  the  effect  of  alloy  scattering  on  material 
properties  and  characteristics  which  determine  device  performance  present 
questions  which  challenge  the  extent  of  this  conclusion. 

The  intent  of  this  paper,  and  other  publications  [7-10],  has  been  to 
offer  a reasonable  calculation  of  the  effects  of  alloy  scattering  on  the 
material  and  transport  properties  of  ternary  and  quaternary  III-V 
semiconductors.  Of  the  three  methods  discussed  for  estimating  the 
magnitude  of  the  scattering  potential,  the  bandgap  difference  is  probably 
least  accurate.  The  other  two  estimations  represent  two  different  views 
of  the  alloy  scattering  potential.  For  some  materials  the  electron 
affinity  difference  technique  and  heteropolar  energy  difference  technique 
give  comparable  magnitudes  for  the  scattering  potential.  For  other  materials 
such  as  Gap,  As  , AlP,  As  and  GaP,  Sb  the  electron  affinity  difference 
gives  a very  small  alloy  scattering  effect  while  the  heteropolar  energy 

i 

difference  gives  a large  alloy  scattering  effect.  Experimental  data  on  these  j 

particular  ternary  alloys  appears  to  be  most  useful  in  experimentally  ; 

determining  which  of  these  models  is  most  accurate  for  the  III-V  semiconductors. 
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Table  1.  Alloy  scattering  potentials  in  electron  volts  for  several 
ternary  III-V  semiconductors  for  a)  the  electronegativity 
theory  b)  the  electron  affinity  difference 

c)  the  energy  band  gap  difference  (AU^^). 


Material 

a)  AUj,jj(x=0.5) 

'>  “eg 

Ga,  In  As 

0.529 

0.830 

1.08 

1-X  X 

InP,  As 

0.581 

0.50 

0.981 

1-X  X 

GaP,  As 

0.637 

0.07 

1.30  ^ 

1-X  X 

InAs,  Sb 

0.801 

0.310 

0.184 

1-X  X 

Al,  In  As 

0.466 

1.320 

1.790 

1-X  X 

AlP,  As 

0.636 

0.08 

0.267 

1-X  X 

Ga,  In  P 

0.559 

0.40 

0.92 

1-X  X 

Al,  In  P 

0.541 

0.90 

1.08 

1-X  X 

Ga_  In  Sb 

0.486 

0.53 

0.515 

1-X  X 

In?,  Sb 

1.32 

0.19 

1.16  5 

1-X  X 

GaP,  Sb 

1.51 

0.06 

1.57 

1-X  X 

I 


r 
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Table  2.  Quaternary  alloy  scattering  parameter  in  electron  volts 

for  quaternary  III-V  semiconductors  for  a)  the  electronegativity 
theory  b)  the  electron  affinity  difference 

c)  the  energy  gap  difference  The  composition  is 

chosen  as  x=y=0.5  for  convenience. 

Material  a)  AU^^CAU^j^)  b)  AU^^iAU^,^)  c)  AUq^CAU^,^) 

0.289  0.263  0.540 

0.536  0.173  0.555 


Ga,  In  P,  As 
1-x  X 1-y  y 


Ga,  In  P,  Sb 
1-x  X 1-y  y 


Al,  In  P,  As 
1-x  X 1-y  y 


0.280 


0.419 


0.581 
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Table  3.  Field-effect  transistor  parameters  based  on  the  model  of  Lehovec 
and  Zuleeg  [25]  for  the  materials  discussed  in  this  paper, 
illustrating  the  effect  of  alloy  scattering. 


Material 

g (mS) 
m 

C (pf) 
gs 

f.j.(GHz) 

r-  (ohms) 
Ds 

f (GHz) 

max 

T.j.(psec) 

GaAs 

30.48 

.3640 

13.33 

445.06 

24.55 

9,25 

Ga  pin  pAs 

a)  No  alloy 

61.39 

.4886 

20.00 

202.31 

35.23 

5.79 

■>>  “en 

54.61 

.4884 

17.79 

228.28 

31.40 

6.51 

48.39 

.4914 

15.67 

247.39 

27.12 

7.39 

43.24 

.4910 

14.02 

278.39 

24.32 

8.27 

a)  No  alloy 

47.55 

.3820 

19.81 

335.42 

39.56 

6.12 

44.16 

.3778 

18.60 

390.17 

38.61 

6.64 

c)  AU_,.,  AU„„ 

35.55 

.3811 

14.85 

456.20 

29.89 

8.25 

a)  No  alloy 

57.85 

.4073 

22.60 

235.35 

41.70 

5.34 

'‘‘'en 

56.73 

.4082 

22.12 

236.15 

40.47 

5.46 

47 . 47 

.4059 

18.62 

294.33 

34.79 

6.49 

42.66 

.4101 

16.56 

303.41 

29.78 

7.28 

Notation:  g = device  transconductance,  C = gate-source  capacitance, 
f^  = gain-bandwidth  product,  small  signal  drain-source 

resistance,  f = maximum  frequency  of  oscillation,  t.j,  = source- 
drain  transit^ime. 

17  -3 

These  calculations  were  nade  for  a device  with  a channel  doping  of  10  cm  and 
the  following  dimensions:  Channel  width  = 0.3ym,  Channel  length  = l.Spm, 
Channel  depth  = 300  gm.  The  gate  voltage  = 0 volts  and  the  drain  voltage 
equals  the  pinch-off  voltage. 


i 
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Figure  1 
Figure  2 

Figure  3 
Figure  4, 

F igure  5 , 

F igure  6 . 


Figure  Captions 

. Electron  scattering  potential  for  ternary  alloys  obtained 

from  the  heteropolar  energy  difference  of  the  binary  constituents. 


Velocity-field  characteristic  for  Ga  In  ^As  obtained  by  Monte 
Carlo  simulation.  The  notation  is:  'Curve  a)  No  alloy  scattering 
p = 9700  cm^/volt 'see;  Curve  b)  Calculated  using  AU  , p=8600 
cm^/volt  sec;  Curve  c)  Calculated  using  AU  , p=86ob  cm^/volt  sec 
Curve  e)  GaAs,  Njj=10^^cm”^ , p=4900  cm^/volt  sec. 

2 

Note:  For  Ga  , „In  ,-As  at  T=300K  p=8450  cm  /volt  sec  [ref. 6]. 

.4/  .00 


Quaternary  alloy  scattering  parameter. 


^(eV^),  for  Ga,  In  P As 
’ 1-x  X 1-y  y 


Velocity-field  characteristic  for  Ga  ^^As  ^ obtained  by 

Monte  Carlo  simulation.  The  notation  is:' 

Curve  a)  No  alloy  scattering;  Curve  b)  Calculated  using 

AU„,  and  AU_.,;  Curve  c)  Calculated  using  AU  . 

LA  EN  Lb 


Velocity-field  characteristic  for  A1  „,-In  „^P  obtained 

• 2^  a/o  a2o  a/D 

by  Monte  Carlo  simulation.  The  notation  is:  Curve  a)  No  alloy 
scattering;  Curve  b)  Calculated  using  AU^.^^;  Curve  c)  Calculated 
using  AUj,^;  Curve  d)  Calculated  using  AU^^. 

Velocity-field  characteristic  for  Ga  m ^b  obtained 

.20./D.OH.lb 

by  Monte  Carlo  simulation.  The  notation  is:  Curve  a)  No  alloy 
scattering;  Curve  b)  Calculated  using  AU^,^;  Curve  c)  Calculated 

using  AUj,j^  and  AU^^. 


APPENDIX  A:  TERNARY  SCATTERING  POTENTIAL  PROM  THE  HETEROPOLAR  CRYSTAL  ENERGY 


III„V 


The  heteropolar  crystal  energy  in  a binary  A B semiconductor  is 


given  in  r.m.k.s.  units  by  [15,16] 


h 

C,„  = T [ — ] exp(-k  R), 

AB  4Tie  ‘•r.  r„  ^ s’ 

o A B 


(A-1) 


where and  are  the  valence  numbers  (3  and  5,  respectively),  r^  and  r^ 
are  the  covalent  radii,  R = 0*5(r^+rg)  is  the  A-B  bond  length,  and  k^  is  the 
Thomas-Fermi  screening  wave  number.  The  factor  b accounts  for  the  fact  that 
the  Thomas-Fermi  approximation  overestimates  screening  for  small  interatomic 
distances  [16].  The  Thomas-Fermi  wave  number  is  given  in  r.m.k.s.  units  by 


, 2 _ 1 .tt,1/3  B 

K - 77  (,— ) 1 

s 4 3 n 

o 


(A-2) 


where  n is  the  valence  electron  density  [15]  and  a„  is  the  Bohr  radius 
o n 


The  valence  electron  density  is  then  given  by 


(A-3) 


where  a is  the  zincblende  lattice  constant, 
o 

If  the  binary  materials  A^^^C^  and  B^^^C  are  alloyed  to  form  C 

and  it  is  assumed  that  the  A-C  and  B-C  bond  lengths  are  equal,  then  the 


fluctuations  in  the  heteropolar  energy  in  the  alloy  can  be  expressed  as  [16] 


- '"ac-'^bcI 


(A-,4) 


.o  bZ  rl 

or  AC  = [ — 

4ne  r, 
o A 


— ]exp(-k  R ) » 
rg  ^ s A 


where  Vegard'r.  law  is  assumed  to  ai>ply  to  both  and  R^ . Thus 


^o  = ^ ’‘^BC 


P-A  = t XTg  + (l-x)r^] 


h 


(A-5) 
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If  this  fluctuation  is  taken  to  be  the  scattering  potential,  then 

AU  = AC,  according  to  the  electronegativity  theory  of  Phillips 
EN 

[15,16].  Equations  1 through  5 are  used  to  calculate  the  entries  in 
Table  1.  Here  the  rationalized  covalent  radii  given  by  Phillips  [23] 
have  been  used  in  the  calculations. 


APPENDIX  B:  QUATERNARY  ALLOY  SCATTERING  RATE 

For  a quaternary  III-V  alloy,  , the  virtual  crystal 

potentials  for  group  III  and  V elements,  respectively,  are  tl3 


"ill  = ’‘"b 


(B-1) 


Uy  = (l-y)U^  + yUp, 

where  U , U_,  U , and  U_  are  the  atomic  crystal  potentials  of  each  element. 
ABC  D 

If  the  A and  B atoms  are  randomly  placed  on  the  group  III  sites  and  the 
B and  C atoms  are  similarly  placed  on  V sites,  the  scattering  rates  due 
to  potential  fluccuations  of  both  III  and  V sites  will  be  proportional 
to  the  square  of  the  matrix  elements  and  the  probability  of  occurrence  of 
each  species  of  atom.  The  probabilities  for  occurrence  in  a random  crystal 
are  (1-x) , x,  (l-y),  and  y for  A,B,C,  and  D atoms,  respectively.  The  matrix 


elements  are  given  by 


|m^l  = /i|j*AU^ij)dr 


(B-2) 


where  AU.  = U^-r  “U.  for  i=A  or  B,  and  AU.-U  -U.  for  i=C,D.  Thus,  the  total 
"illli  ivi 


scattering  rate  for  the  quaternary  alloy  is  given  by 

— -x-  xU-x)  + y(i-y)|My|^ 

QA 


(B-3) 


l^iyl^  = 1/<,*(U^-Up)i|;drl' 


(B-‘+) 


25 


If  U.-U  is  regarded  as  the  change  in  group  III  potential  with  group  V 
A D 

atom  fixed,  and  is  regarded  as  the  change  in  group  V potential  with 

group  III  atom  fixed,  then  the  composition  weighted  averages  for 


U^-Ub  and  are 


Ub-Ua  = (l-y)(U„^-U,^)  t y(u„„-u^^) 


BC  AC 


BD  AD' 


(B-5) 


U^-Ub  = (i-xXUbb-Ubc)  - -(Uab-Uac) 


Here  (U  -U.„)  represents  the  difference  in  the  effective  potential  for 
dC  AC 

electrons  with  either  group  B or  group  A atoms  on  group  III  sites  with  a 
C atom  definitely  on  a group  V site,  and  so  on  for  the  other  three  terms. 

By  substitution  of  Eqn,  (5)  into  Eqn.  (4)  and  neglecting  the 
"overlap"  integrals  one  obtains 

QA 
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Here  the  AU's  on  the  r.h.s.  of  Eqn.  7 are  ternary  scattering  potentials. 

For  example,  AU.  is  the  scattering  potential  of  the  ternary  A B D 

and  is  the  scattering  potential  of  the  ternary  y^y'  addition, 

each  material  parameter  in  Eqn.  7 is  a material  parameter  of  the  quaternary 
alloy,  which  can  be  estimated  from  binary  and  ternary  material  parameters 
by  the  interpolation  procedures  previously  described  [10]. 


o>  CO 

o o o o o o 


(A9)*^^nV  ‘lVllN310d  ONIdBllVOS  AOHIV 


0.2  0.4  0.6  0.8 

ALLOY  COMPOSITION  , X 


=10  cm  T=300K 


ELECTRIC  FIELD  INTENSITY  (kV/cm) 


IAUqaI  , Alloy  Scattering  Parameter  (eV 
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rvj 


In  As 


ELECTRIC  FIELD  INTENSITY  ( kv/cm) 


ELECTRIC  FIELD  INTENSITY  (kV/cm) 


ELECTRIC  FIELD  INTENSITY  (kv/cm) 
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NEGATIVE  RESISTANCE  AND  PEAK  VELOCITY 
IN  THE  CENTRAL  (000)  VALLEY  OF  III-V  SEMICONDUCTORS* 

J.  R.  Hauser,  T.  H.  Glisson  and  M.  A.  Littlejohn 
Electrical  Engineering  Department 
N.  C.  State  University 
Raleigh,  N.  C.  27607 

ABSTRACT 

The  negative  resistance  in  III-V  materials  such  as  GaAs  at  large 
electric  fields  is  generally  recognized  as  arising  from  the  transfer 
of  electrons  from  the  central  (000)  valley  to  higher  lying  minima  in 
the  conduction  band.  Monte  Carlo  transport  studies  show  that  the 
negative  resistance  effect  is  still  present  in  III-V  materials  when  the 
valley  spacing  is  increased  to  large  values  (>  0.5  eV)  and  even 
present  when  the  higher  minima  are  eliminated  entirely  from  the 
calculations.  This  negative  resistance  arises  from  basic  transport 
properties  of  the  central  valley.  Studies  are  presented  of  the  basic 
negative  resistance  effect  in  the  central  valley  of  III-V  materials  as 
well  as  studies  of  Al,  In  As(x'''0.75)  and  Ga,  In  As  (x'^0.6)  which  are 
two  specific  materials  where  the  negative  resistance  effect  is  due 
predominantly  to  the  central  valley. 


■•'•This  worr  was  supported  by  a research  grant  from  the  Office  of  Naval  Research, 
Arlington,  VA. 
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NEGATIVE  RESISTANCE  AND  PEAK  VELOCITY 
IN  THE  CENTRAL  (000)  VALLEY  OF  III-V  SEMICONDUCTORS 

I.  Introduction 

The  negative  resistance  effect  as  first  reported  by  Gunn  [1,2] 
has  been  observed  in  most  of  the  III-V  binary  and  ternary  compounds 
which  are  direct  bandgap  materials.  The  major  factors  determining 
the  negative  resistance  characteristic  are  now  thought  to  be  fairly 
well  understood,  and  a good  summary  has  recently  been  given  by  Ridley  [3]. 

A brief  summary  of  the  conventional  understanding  is  useful  in  relationship 
to  the  present  work. 

The  negative  resistance  effect  in  GaAs  has  been  identified  with  a 
transfer  of  electrons  at  large  electric  fields  from  the  high  mobility 
r or  central  valley  to  higher  lying  L or  X minima  having  lower  mobilities. 
The  validity  of  this  basic  model  has  been  verified  by  simple  analytical 
calculations  [3-5]  as  well  as  by  more  detailed  Monte  Carlo  calculations 
[6-11].  The  energy  bands  for  a general  III-V  semiconductor  are  shown  in 
Figure  1,  where  is  the  energy  separation  between  the  central  valley  and 
the  next  lowest  valley,  which  is  assumed  to  be  the  L valley.  In  some 
materials  a three  band  model  including  both  the  X and  L minima  is  required 
to  give  good  agreement  between  theory  and  experimental  results. 

In  tr.e  presence  of  a large  electric  field,  central  valley  electrons 
are  heated  by  the  field  to  large  kinetic  energies,  and  as  the  average 
kinetic  energy  approaches  the  valley  separation,  a large  percentage  of  the 
electrons  transfer  to  the  upper  valley.  This  transfer  is  enhanced  by  the 
larger  effective  mass  of  the  upper  valleys  as  compared  with  the  central 
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valley.  Ridley  has  used  an  energy  balance  approach  to  estimate  the  average  ] 

energy  of  the  electrons  as  a function  of  electric  field  and  equated  this 
average  energy  to  the  valley  separation  to  estimate  the  threshold  electric  I 

field  for  a variety  of  transferred  electron  materials  [3].  | 

Simple  considerations  of  the  transferred  electron  effect  lead  to  the  | 

conclusion  that  high  peak  velocities  should  be  achieved  in  materials  with  j 

a)  a large  low  field  mobility  and  b)  a large  energy  separation  between  | 

I 

1 

the  central  valley  and  the  higher  lying  minima.  The  binary  semiconductor  i 

GaAs  has  a relatively  small  valley  separation  on  the  order  of  0.31  eV.  j 

The  separation  for  InP  is  somewhat  larger  at  0.60  eV  and  InP  is  predicted  | 

to  have  a larger  peak  velocity  than  GaAs  although  the  field  at  which  the  j 

peak  velocity  occurs  is  alsc  somewhat  larger. 

The  largest  valley  separation  in  the  binary  III-V  materials  is  1.11  eV 

! 

in  InAs.  This  large  valley  separation,  however,  cannot  be  used  in  TED  1 

devices  because  of  the  low  breakdown  voltages  of  InAs  due  to  the  small 

bandgap  of  0.36  eV,  Thus,  of  the  binary  III-V  materials,  InP  is  predicted 

to  have  the  largest  peak  velocity  of  the  materials  which  are  useful  for  ; 

TEDS. 

The  ternary  and  quaternary  III-V  materials  provide  a wider  range  of 
bandgaps  and  valley  separation  for  potentially  higher  peak  velocities 
than  can  be  achieved  in  the  binary  materials.  The  quaternary  Ga^  y^^y 

lattice  matched  to  InP  substrates  has  recently  been  predicted  to  have  a 
peak  velocity  larger  than  that  of  either  GaAs  or  InP  [12].  The  large  low 
field  mobility  needed  for  high  velocities  requires  low  bandgap  materials 
with  small  effective  nasses.  However,  too  low  a bandgap,  as  in  InAs,  leads  : 

to  low  breakdown  voltages.  Thus  a compromise  must  be  found  between  these 
two  requirements.  From  breakdown  voltage  and  device  considerations,  the 
bandgap  of  high  velocity  materials  for  either  TEDs  or  FETs  should  probably 


be  around  1 eV  or  larger. 
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If  a mininium  bandgap  of  around  1 eV  is  selected  then  all  the  III-V 

ternary  and  quaternary  materials  can  be  searched  to  find  which  material 

has  the  largest  valley  separation.  This  material  turns  out  to  be  Al^_^In^As 

for  which  the  energy  band  diagram  vs.  composition  is  shown  in  Figure  2. 

At  a composition  of  x=0.75  the  energy  gap  is  about  0.91  eV  while  the  F-to-L 

valley  separation  is  about  1.12  eV.  Monte  Carlo  calculations  for  this 

composition  (to  be  discussed  in  the  next  section)  lead  to  a low  field 

2 

mobility  of  about  11,700  cm  /V.sec  and  this  coupled  with  the  large  valley 
separation  leads  to  an  expected  large  pea)<  velocity  before  the  onset  of  a 
negative  resistance  effect. 

Another  ternary  material  which  has  a large  valley  separation  and  also 
has  a large  bandgap  is  Ga^  ^In^As  with  x'^^O.  3-0.4.  For  example  at  x=0.4 
the  bandgap  is  0.86  eV  and  the  valley  separation  is  0.72  eV.  These  values 
are  not  quite  as  good  as  those  for  the  Al^  ^In^As  system  but  the  technology 
for  Ga^_^In^As  is  much  more  highly  developed.  The  quaternary  Ga^  x^'^x^l  y^^y 
previously  reported  on  [12j  also  has  a favorable  bandgap  and  valley  separation 
for  large  peak  velocities. 

Monte  Carlo  calculations  on  both  the  Al,  In  As  and  Ga,  In  As  systems 

1-X  X 1-X  X 

have  not  shown  as  large  a pe-k  velocity  as  initially  expected.  This  has 
been  found  to  be  due  to  fundamental  physical  limitations  for  the  peak 
velocity  which  occur  within  the  central  valley.  For  valley  separations 
larger  than  about  0.5  eV , the  peak  velocity  and  the  threshold  field  have 
been  found  to  be  determined  almost  entirely  by  the  properties  of  the 
central  valley.  Materials  with  large  valley  separations  will  be  referred 
to  in  this  work  as  central  valley  dominated  materials  since  such  materials 
show  a peak  velocity  and  a negative  resistance  determined  mainly  by  the 
central  valley.  The  high  field  properties  of  such  materials  are  discussed 


in  detail  in  the  next  section. 


Figure  3.  The  material  parameters  used  in  the  calculations  were 
obtained  from  the  binary  III-V  parameters  as  discussed  in  previous 
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publications  [10-12]  and  are  given  in  Table  I.  Included  in  the 

calculations  are  scattering  processes  due  to  acoustic  phonon,  polar 

16  3 

optical  phonons,  ionized  impurity  (10  /cm  impurity  density),  piezo- 
electric, alloy  scattering  processes,  and  equivalent  and  nonequivalent 

7 

intervalley  processes.  The  calculated  peak  velocity  is  about  2.7x10  cm/sec 
at  a field  of  about  4000  V/cm.  This  peak  value  is  about  35%  higher  than 
the  calculated  and  measured  values  for  GaAs  and  this  verifies  to  some 
extent  the  discussion  of  the  previous  section  as  to  the  need  for  a 
large  valley  separation  to  achieve  large  peak  velocities.  However, 
the  peak  velocity  is  not  as  large  as  was  expected  before  the  calculations 
were  performed. 

The  numerical  values  shown  at  various  points  along  the  curve  of 

Figure  3 give  the  percentages  of  electrons  in  the  L plus  X minima.  At 

the  peak  of  the  velocity  curve  it  is  seen  that  only  0.07%  of  the  electrons 

4 

have  been  transferred  to  the  upper  valleys.  At  10  V/cm  which  is  f.ar  into 
the  negative  resistance  region  only  9.2%  of  the  electrons  have  been 
excited  to  the  upper  valley.  These  percentages  which  were  obtained  from 
the  Monte  Carlo  calculations  are  much  too  low  to  account  for  the  peak  in 
velocity  and  the  negative  resistance  for  this  material.  For  example  at 

4 

10  V/cm  a transicr  of  all  electrons  back  to  the  central  valley  would 
increase  the  velocity  by  no  more  than  10%  and  this  is  much  too  small  an 
effect  to  eliminate  the  negative  resistance.  This  leads  to  the  conclusion 
that  the  peak  velocity  is  being  controlled  by  the  central  valley  in  this 


material  and  is  not  due  to  electror,  transfer  to  the  upper  valley. 


The  dominance  of  the  central  valley  in  determining  the  peak  velocity 
is  shown  in  Figure  4 which  compares  the  calculated  velocity-field 
relationship  of  Figure  3 with  that  obtained  by  including  only  the 
central  valley  in  the  Monte  Carlo  calculations.  The  two  calculations 
give  very  similar  velocity-field  relationships  with  the  peak  velocity 
being  unchanged  by  the  elimination  of  scattering  to  the  upper  valleys. 

The  solid  curve  gives  an  inherent  negative  resistance  phenomena  which 
has  not  previously  been  recognized  to  occur  for  III-V  materials  within 
the  central  valley  alone.  A comparison  of  the  two  curves  in  Figure  4 
shows  that  transfers  to  the  upper  valleys  increase  the  magnitude  of  the 
negative  resistance  beyond  the  peak  but  the  upper  valleys  are  not  required 
for  the  existence  of  the  negative  resistance.  Upper  valley  transfers 
are  also  seen  to  cause  a larger  drift  velocity  at  very  large  fields  (at 
10^  V/cm  for  example). 

The  dominance  of  the  central  valley  is  not  limited  to  just  A1  In  As 

but  has  also  been  seen  in  other  ternary  materials  such  as  Ga  In  As  with 

; X 0.4  - 0.6.  The  velocity  field  curve  for  this  ternary  with  x=0.6  is 

j shown  in  Figure  5.  This  particular  material  has  an  energy  gap  of  0.65  eV 

I and  a valley  separation  of  0.90  eV . As  with  the  A1  25^^  75^^  case,  the 

! 

presence  of  the  upper  valleys  has  essentially  no  effect  on  the  peak 
velocity  or  threshold  field  for  negative  resistance.  The  upper  valleys 
again  increase  the  negative  resistance  effect  and  increase  the  high  field 
velocity . 

The  two  materials  discussed  so  far  have  large  valley  separations, 


and  this  is  the  condition  for  the  peak  velocity  and  threshold  field 
being  determined  by  the  central  valley.  For  lower  valley  separations, 
however,  such  as  the  0.31  eV  for  GaAs  the  question  arises  as  to  how 
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important  is  the  intervalley  transfer  of  electrons  in  determining  the 
peak  velocity.  Calculations  for  GaAs  have  shown  that  without  the  L and 

7 

X valleys  present  the  velocity  would  peak  at  about  2.5  x 10  cm/sec  as  opposed 

7 

to  the  value  of  about  2.0  x 10  cm/sec  including  the  upper  valleys.  Thus 
the  transfer  of  carriers  to  the  upper  valleys  reduces  the  peak  velocity  in 
GaAs  by  about  20%.  A somewhat  similar  central  valley  dominance  of  the  velocity 
field  relationship  has  been  seen  in  Monte  Carlo  calculations  on  CdTe  [13]. 

This  material  also  has  a large  valley  separation  similar  to  that  of  the 
III-V  materials  discussed  above. 

The  central  valley  dominated  materials  such  as  shown  in  Figure  4 are 
predicted  by  the  Monte  Carlo  calculations  to  have  a large  peak-to-valley 
ratio  for  the  velocity-field  curve.  The  peak-to-valley  ratio  in  Figure  4 
ranges  from  about  4 to  7 for  peak  fields  from  20  kV  to  100  kV.  For  GaAs 
with  a smaller  valley  spacing,  the  peak-to-valley  ratio  is  at  most  about 
2.5  by  both  experimental  measurement  [14]  and  theoretical  calculations 
[15].  This  large  peak-to-valley  ratio  for  the  central  valley  dominated 
semiconductors  may  be  especially  useful  in  certain  device  applications. 

The  preceding  discussion  has  demonstrated  that  when  the  valley  spacing 
becomes  large  the  peak  velocity  and  threshold  field  are  no  longer  determined 
by  electron  transfer  from  the  central  valley  to  higher  lying  minima.  The 
pea):  velocity  is  rather  determined  by  the  fundamental  transport  properties 
of  the  central  valley.  In  order  to  understand  these  central  valley 
limitations  a study  has  been  made  of  the  high  field  properties  of  model 
semiconductors  with  only  a central  valley.  These  studies  are  discussed 


ir  detail  in  the  next  section. 
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III.  Transport  in  Central  Valley  Only 


r- 

. 

This  section  considers  a semiconductor  in  which  the  spacing  between 
the  central  valley  [iX 000)  valley]  and  higher  lying  minima  (either  X or  L 
valley)  is  so  large  that  only  the  central  valley  must  be  considered  in 
the  transport.  In  all  of  the  III-V  ternary  and  quaternary  materials 
the  spacing  can  never  be  made  so  large  that  it  has  no  effect  on  the 
trcinsport  process.  However,  as  the  previous  section  illustrates  the 
spacing  can  be  so  large  that  intervalley  transfers  have  a negligible 

I effect  on  the  peak  velocity  and  threshold  field  for  certain  materials. 

I 

[ A study  of  transport  in  only  the  central  valley  can  thus  lead  to  an 

> 

i 

I understanding  of  the  physical  processes  leading  to  a peak  velocity  and 

I negative  resistance  effect  in  these  materials  as  well  as  provide  an  upper 

1 

I limit  to  the  peak  velocity  for  other  materials.  Since  the  A1  .,j-In  ..j-As 

I • ^ D • / D 

comes  the  closest  of  all  the  ternary  III-V  materials  to  being  dominated 

1 

[ by  the  central  valley,  the  parameters  for  the  calculations  presented  in 

this  section  using  only  the  central  valley  have  been  selected  as  those 
of  this  material  which  are  listed  in  Table  I , 

The  negative  resistance  and  peak  velocity  in  the  central  valley 
have  been  founc  to  be  due  to  the  fundamental  properties  of  polar  optical 
scattering  in  this  valley.  This  is  shown  in  Figure  6 which  shows  the 
calculated  velocity-field  curve  (solid  curve)  for  only  polar  optical 
scattering  in  the  central  valley.  The  dotted  curve  shows  the  velocity- 
field  curve  when  all  central  valley  scattering  processes  are  included. 

Polar  optical  scattering  alone  is  seen  to  determine  the  general  shape  of 
the  curves  with  the  other  scattering  processes  simply  reducing  the  velocity 
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optical  ;;ca  I Icri  is  to  he  expelled  ;,iiii;e  it  is  the  luijor  eiier'f.y  lo.ss 
process . 

The  central  valley  negative  resistance  for  cinly  pol.ir  optical  scattering 

can  be  understood  with  fairly  simple  models.  In  steady  state  the  use  of 

energy  balance  and  jmentum  b.ilancc  exprera. i ons  gives 

rift  = m'-'-v , (1) 

m 

A£(P  -P  ) 

oEv  = — , (2) 

T 

C 

where  m"  is  the  effective  mss,  v is  the  average  drift  velocity,  is  the 
rromentum  relaxation  time,  is  the  energy  relaxation  time,  AS  = is  the 

energy  lost  or  gained  in  .a  polar  optical  scattering  event,  is  the  probability 
of  phonon  er.ission  anci  ! is  the  pro'  oility  of  phonon  absorption.  At  large 

<3 

energies 


r « ^ 

‘ e * a 


-)  = tan:.(x/2). 


where  >:  = d'lOii^^/kT . Solving  both  Equat'ens  (1)  ind  (2)  for  the  velocity  gives 


d'n...  tanhidl^-  /2KT)  , 

to  to  i 


If  is  independent  of  field  then  Equation  (4)  ju'eiicts  a linearly  increasing 
velocity  v.-ith  fiel,:  ..ucn  as  is  observed  in  figure  6 in  the  positive  lesistar.re 


region. 


i:.  Equation  (5)  is  constant  or  varies  slowly  with  field  then 


this  equa:  n predicts  a decreasing  velocity  with  field  such  as  observed  in 
figure  f in  the  negative  resistance  region. 
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The  agreement  between  these  two  expressions  and  the  Monte  Carlo 
calculations  is  shown  in  Figure  7.  The  solid  curve  is  the  Monte  Carlo 
results  while  the  dotted  curves  show  Equations  (4)  and  (5)  in  the  low 
field  and  high  field  regions.  In  using  Equations  (4)  and  (5)  and 
have  not  been  taken  as  constant  but  the  curve  of  t si  ^wn  in  Figure  7 has 

been  used  in  the  calculations.  This  t was  obtained  from  the  Monte  Carlo 

. . a 

calculations  and  is  the  mean  time  between  polar  optical  scatterings  . 

Simple  analytical  expressions  for  or  as  a function  of  field  are 
not  easily  derived  so  the  results  of  the  Monte  Carlo  calculations  were 
used  to  verify  Equations  (4)  and  (5). 

The  good  agreement  of  the  Monte  Carlo  calculations  with  Equations  (4) 
and  (5)  in  the  limiting  regions  shows  that  at  low  fields  the  velocity  is 
primarily  determined  by  momentum  balance  considerations  while  at  high 
fields  in  the  negative  resistance  region  the  velocity  is  determined  Dy 
energ"  balance  considerations.  With  optical  phonon  processes  there  is 
a fixed  energy  loss  per  scattering  event  and  the  drift  velocity  must 
decrease  with  field  or  too  much  ei.ei'gy  will  be  gained  from  the  field  to 
be  lost  in  the  scattering  events. 

An  estination  of  the  peak  velocity  which  can  be  achieved  by  carriers 
in  the  central  valley  can  be  obtained  by  equating  the  limiting,  expressions 
as  given  by  Equations  (4)  and  (b)  and  solving  for  qEi/m*.  This  gives 
for  the  peak  velocity 

a.  The  use  of  the  same  t for  both  i and  t will  be  approximately  valid 
as  long  as  only  polar  optical  scattering  is  important.  When  other 
scattering  events  are  present  and  will  be  significantly  different. 


Values  of  v for  several  binary  and  compound  materials  are  shown  in 
max 


Figure  8.  Also  shown  for  several  materials  are  the  peak  velocity  values 
which  result  from  Monte  Carlo  calculations  including  all  bands  and  all 
scattering  processes.  For  the  binary  materials  GaAs , InP  and  InAs  the 
Monte  Carlo  results  give  peak  velocities  which  are  about  82%-go%  of  the 
values  predicted  by  Equation  (6).  For  the  two  ternary  and  one  quaternary 
materials  shown,  the  Monte  Carlo  results  differ  from  Equation  (6)  by  large 
amounts.  The  lower  Monte  Carlo  results  for  the  ternary  and  quaternary 
materials  is  due  to  the  presence  of  alloy  scattering  in  these  materials 
which  is  not  present  in  the  binary  materials.  If  alloy  scattering  had 
not  been  included  all  of  the  Monte  Carlo  results  would  have  been  around 
80%-90%  of  the  limiting  value  given  by  Equation  (6). 

Monte  Carlo  results  are  not  shown  for  InSb  and  GaSb . The  bandgap 
of  InSb  is  too  low  for  this  material  to  be  used  near  the  peak  velocity. 

For  GaSb  the  presence  of  the  L mininv-i  at  only  0.09  eV  above  the  F minima 
causes  the  negative  resistance  effect  not  to  occur  in  this  material.  The 
data  in  Figure  6 indicates  that  there  is  a general  trend  toward  larger 
peak  velocities  as  the  bandgap  of  the  semiconductor  decreases. 

In  addition  to  the  fundamental  role  of  polar  optical  scattering  in 
the  central  valley,  the  band  nonparal^olicity  also  has  an  influence  on 
the  peak  velocity  which  can  be  achieved.  Figure  9 shows  the  velocity  field 
curves  for  the  central  valley  only  of  A1  yt,As  with  different  a values 


wher( 
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eci  t C.6)  = Y(e). 


(7) 


The  theoretical  value  of  a is 


1 ,,  m*.2 

a = ^ ( 1-  — ) , 

o 


(8) 


,-l 


and  for  A1  25^^  75^^  value  is  1,04  eV  corresponding  to  one  of  the 
curves  in  Figure  9.  Also  shown  are  velocity  field  curves  for  a values 
an  order  of  magnitude  smaller  and  an  order  of  nagnitude  larger  than  the 
theoretical  value.  For  values  of  a less  than  0,1  the  calculated  values  are 
close  to  the  0.104  curve  of  Figure  9.  At  the  theoretical  value  of  a 
the  nonparabolicity  is  seen  to  reduce  the  low  field  mobility  by  about 


37%,  to  reduce  the  peak  velocity  slightly,  and  to  increase  the  tlireshold 
field  from  around  3 KV/cm  to  around  4 KV/cm.  A very  large  nonparabolicity 
(a=10.4  eV  ^) , which  fortunately  does  not  occur  in  the  III-V  materials, 
is  seen  to  very  drastically  reduce  the  low  field  mobility  and  the  peak 
velocity.  The  nonparabolocity  will  be  most  important  in  materials  with 
narrow  bandgaps  and  large  valley  spacings,  since  a is  large  for  such 
materials  and  the  carriers  become  heated  to  large  energies  before  inter- 
valley transfer  occurs. 
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IV.  Summary  and  Conclusions 


The  presently  accepted  theory  of  negative  resistance  effects  in 
III-V  semiconductors  attributes  the  effect  to  the  transfer  of  electrons 
from  a high  mobility  central  ( D valley  to  lower  mobility,  higher  lying 
energy  valleys  (L  or  X),  The  present  work  has  shown  that  in  materials 
with  a large  valley  spacing  the  negative  resistance  becomes  dominated 
by  the  central  valley.  Monte  Carlo  calculations  have  shown  that  the 
presence  of  the  upper  valleys  is  not  required  for  a negative  resistance 
effect  but  that  polar  optical  scattering  acting  in  the  central  valley 
alone  gives  rise  to  a peak  velocity  and  negative  resistance  effect.  For 
materials  with  valley  spacings  of  0.5  eV  or  larger  the  peak  velocity  and 
threshold  field  has  been  found  to  be  determined  almost  entirely  by  the 
central  valley  negative  resistance  effect.  Transfers  to  the  upper  valleys 
influence  mainly  the  magnitude  of  the  negative  resistance  and  the  carrier 
velocity  near  the  valley  minimum  of  the  velocity  field  curve. 

Materials  which  have  large  valley  spacings  and  are  dominated  by  the 

central  valley  have  large  peak-to-valley  velocity  ratios  as  predicted 

by  Monte  Carlo  calculations.  The  ternary  materials  Al,  In  As  (x^-0.75), 

^ 1-x  X 

Ga,  In  As(x'v<0.5)  and  the  quaternary  Ga,  In  P,  As  (x^-.VS,  y'v.4)  are 
1— X X ^ 1“X  X I'-y  y 

typical  of  central  valley  dominated  materials.  The  largest  peak  velocities 
which  have  so  far  been  calculated  by  the  Monte  Carlo  technique  have  been 
found  in  these  materials. 
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.^gAs  Material  Parameters 

Used  in  the  Calculations 

A.  Bulk  Material 

Parameters 

Parameter 

Value 

Parameter 

Value 

Lattice  Constant 

5.959xl0“®cm 

LO  Phonon  energy 

0.03605  eV 

Density 

3 

5.27  gm/cm 

Sound  velocity 

5 .llxlO^cm/sec 

Electron  Affinity 

Diff.  1.32  eV 

Optical  Dielectric 
Constant 

10.78 

Piezoelectric  Constant  .056  coulomb /m 

Static  Dielectric 
Constant 

13.46 

B.  Valley-Dependent  Material  Parameters 


Conduction  Band  Valley 


Parameter 

rxooo) 

X(IOO) 

L(lll) 

Acoustic  Deformation  Potential  (eV) 

6.73 

9.35 

8.83 

Effective  Mass  (in*/in  ) 

0 

0.028 

0.62 

0.34 

Non-Parabolicity  (eV 

1.04 

0.204 

0.814 

Energy  Band  Gap  (eV) 

0.909 

2.113 

2.025 

(relative  to  valence  band) 

Optical  Deformation  Potential  (eV/cm) 

0 

0 

3x10® 

Optical  Phonon  Energy  (eV) 

- 

- 

0.0336 

Intervalley  Deformation  Potential 
(eV/cm) 

From  r(OOO) 

0 

6.25x10® 

5x10® 

From  X(IOO) 

6.25x10® 

6.25x10® 

5x10® 

F rom  L( 111) 

5x10® 

5x10® 

5x10® 

Intervalley  Phonon  Energy  (eV) 

From  r(OOO) 

0 

0.0253 

0.0262 

From  X(IOO) 

0.0253 

0.0253 

0.0309 

Fro-  1(111) 

0.0262 

0.0309 

0.0254 

Numl;er  of  Equivalent  Valleys 

1 

3 

u 
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Figure  Captions 


Figure  1.  Energy  band  diagram  lor  general  III-V  direct  bandgap  semiconductor. 
Figure  2.  Calculated  energy  gaps  vs.  composition  for  Al^_^ln^As. 

Figure  3.  Calculated  velocity-field  curve  for  A1  25^^  75^^'  values  along 

the  curve  show  the  percentages  of  electrons  in  the  upper  ( 111  and 
100)  valleys. 

Figure  4.  Comparison  of  velocity-field  curves  for  A1  25^^  75^^  using  central 
valley  only  and  using  all  valleys. 

Figure  5.  Velocity-field  curves  for  Ga  ^^In  gAs. 

Figure  6.  Velocity-field  curve  for  central  valley  with  just  polar  optical 
scattering . 

Figure  7.  Comparison  of  Monte  Carlo  calculations  with  simple  momentum  and 
energy  balance  expressions. 

Figure  8.  Calculated  upper  limits  to  peak  velocity  considering  central 
valley  only. 

Figure  9.  Velocity  Field  curves  for  A1  25^^  75^^  with  different  nonparabolicity 
factors.  The  theoretical  value  is  a = 1.04  eV 
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Velocity-field  characteristics  of  Gai_;,lnj,Pi_yASy  quaternary 
alloys* 

M.  A.  Littlejohn,  J.  R.  Hauser,  and  T.  H.  Glisson 

Department  of  £/ectrlcal  Engineering,  .Vurth  Caroiina  Stote  L'niverMfy.  Halcigh.  \orih  Curulinu  27607 
(Received  18  October  1976;  in  final  form  27  December  1976) 

The  electron  drift-veliKity-electric-field  relationship  has  been  calculated  for  the  Ga,  ,ln,l*  As, 
quaternary  alloy  using  the  Monte  Carlo  method  Emphasis  has  been  placed  on  the  compositional  range  for 
which  the  alloy  is  lattice  matched  to  GaAs  and  InP  These  calculations  suggest  that  this  qualeriiary  offers 
promise  as  a material  for  microwave  semiconductor  devices,  including  field-elfect  transistors  and 
transferred  electron  devices. 


PACS  numbers.  72  20.Er,  8.‘i  .lO.De,  72  20  Dp,  8.8..)0  Tv 

Recent  interest  in  the  Ga,.,In,Pj.^.ASy  quaternary  alloy 
has  centered  around  the  applications  of  this  material 
for  optical  devices,  such  as  photoemission  cathodes, 
double-heterojunction  lasers,  and  litiht-emitting 
diodes. One  major  advantage  of  the  quaternary  alloy 
system  for  these  applications  is  that  the  material  can 
be  synthesized  with  a fixed  lattice  constant  which  is 
matched  to  a substrate  over  a wide  enerjry  band-tjap 
range."  This  fact  is  also  important  in  material  consid- 
erations for  microwave  devices,  such  as  the  metal- 
epitaxial- semiconductor— field -effect  transistor 
(MFSFET)  and  transferred  electron  devices.  The  pres- 
ent study  has  been  concerned  with  the  electron  drift- 
velocity— electric-field  relationship  in  Ga,.,In,P,.yAs, 
quaternary  alloys  which  can  be  lattice  matched  to  either 
GaAs  or  InP.  The  results  show,  especially  for  lattice 
matching  of  the  quaternary  to  InP,  that  substantial  im- 
provements in  low  field  drift  mobility  and  peak  drift 
velocity  can  be  obtained  in  comparison  to  GaAs,  InP, 
and  the  Ga,In,.,A.s  ternary  system,  which  are  all  im- 
piirtant  materials  for  microwave  devices. 

The  computer  simulation  of  transport  processes  used 
in  obtaining  the  results  presented  here  has  been  de- 
scribed in  previous  publication;?’-''  and  is  similar  to 
other  such  programs  discussed  in  the  literature.'’"" 

The  allowed  scattering  mechanisms  in  the  transport 
simulation  include  acoustic  phonon  scattering,  optical 
phonon  scattering,  piezoelectric  scattering,  equivalent 
and  nonequivalent  intervalley  scattering,  ionized  im- 
purity scattering,  and  random  potential  alloy  scatter- 
ing." The  program  allows  for  calculations  with  any  or 
all  of  these  mechanisms  in  either  the  rfOOO).  .Y(IOO), 
or  Edll)  nonparabolic  conduction  bands. 

Applied  Pbysics  Letters,  Vol.  30,  No  5,  1 March  1977 


The  manner  in  which  the  material  parameters  for  the 
quaternary  alloy  are  calculated  as  a function  of  alloy 
composition  (v,  v)  is  as  follows.  First,  the  accepted 
values  of  binary  material  parameters  for  GaAs,  ' InP,"-'* 
InAs,"-'-'  and  GaP  arc  used  in  theoretical 
models'"-"  to  calculate  the  material  parameters  tor  the 
four  possible  ternary  combinations  of  tliese  (our  binary 
materials.  For  some  parameters,  such  as  the  equiva- 
lent and  nonequivalent  intervalley  deformation  poten- 
tials, such  theoretical  models  do  not  exist,  and  in  these 


,,  COMPOSITION  M ONO  I al  TICE  - M6TCHF  0 LINES 


I'KI.  1.  Ixiw-fu'UI  mobility  of  ('.:i|_,^n,l’|^ASj,  l.-itlii-i-  ni  iti-lu-d 
to  Ini’  1111(1  (i.-iAs.  Till- (|ii.-it('i-ii.ii-y  (lo|iiiig  level  is  111''  (-nT  . 
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(8)  6614 
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Eg(EV)  E^JCV) 

loe  016 

103  850 

972  661 


Go,.,ln,P,^  *•, 

Lo^iic* 'molchvd  to  InP 


ELECTRIC  FIELD  INTENSITY  (kV/cm) 


FIG.  2.  Vclocity-fielil  cli:iracteristic  »[  C'i:ii.^ii,l’|^A.Sj,  liiltici 
matched  to  Ini’.  The  dopint;  level  i.s  In*''  cm"'. 


cases  linear  interpolations  are  used  for  the  ternary 


parameters."'"  From  the  ternary  material  parameters, 


the  following  equation  is  used  to  calculate  a given  qua- 
ternary material  parameter,  Q,  as  a function  of  the 
quaternary  alloy  composition  (v,  v): 


Q(^,  v)=  [rd  - .v)[(l  - v)r,.,{.v)  + yT„(.v)|  ^ v(l  - v) 

X[(l  - v)r,,(v)  - vT,3(  v)|}lv(l  - v)4  v(l  - v)|-*, 


where  T, ^ is  the  material  parameter  (such  as  band  gap, 
deformation  potential,  etc.)  for  the  ternary  alloy  of 
binary  materials  i and  i.  In  this  paper,  GaP  is  material 
1,  InP  is  material  2,  InAs  material  3.  and  GaAs  mate- 
rial 4.  Thus  7,2(v)  is  a material  parameter  for 
Gaj.,In,P  Equation  (1)  is  an  empirical  relation  which 
reduces  to  the  correct  ternary  expression  (i.e. , for  v 
or  A'  equal  to  either  0 or  1),  and  in  the  limit  reduces  to 
the  correct  expression  for  the  binary  materials.  Also, 
for  V - v = 0.  5 Eq.  (I)  reduces  to  the  average  of  the 
four  ternarv  alloy  parameters.  For  other  \ and  v 
values  a smooth  curve  is  generated  which  correctly 
matches  tlie  ternarv  boundary  parameters.  This  empir- 
ical procedure  for  choosing  the  quaternary  material 
parameters  is  similar  to  that  used  by  Moon  cl  for 
lattice  constant  and  energy  band  gap. 


Microwave  device  technology,  such  as  that  for  the 
MKSFET,  usually  requires  thin  epitaxial  lavers  grown 
on  a suitable  substrate.  The  work  discussed  here  is 


1 1 ] 


Wdhoal  0<>0y  SCOtlC'io^ 


O'lOy 


/ 


FLiCS'C  riElO  intensity  {«v  cml 


I lC'.  a.  \'ck)vily-lielil  clKiniclcri.slic  "f 

wifii  iind  without  riindoni  potcntiiil  alloy  scattering.  Shown  for 
coTniiarison  are  die  velocity-riehl  curves  lor  G.aY\s,  Iid’,  and 
Gaii.f.Iiio^jAs.  Tile  doping  level  is  ID*'  ein"*. 


limited  to  the  quaternary  alloys  which  are  lattice 
matched  to  either  GaAs  or  InP,  since  these  materials 
are  presently  the  most  popular  binary  Etl-V  materials 
used  for  substrates  in  microwave  applications.  The 
values  of  v and  r required  for  lattice  matching  to  InP 
(lattice  constant  - 5.809  A)  and  GaAs  (lattice  constant 
' 5.042  Al  can  be  found  in  Fig.  1 of  Ref.  6.  This  figure 
has  been  reproduced  in  this  work  using  Eq.  (1).  These 
values  of  v and  \'  for  lattice  matching  are  given  to  a 
good  approximation  by  the  linear  equations 


v = 2.  197(1  - v),  matching  to  InP 


v = l -2.  179i 


matching  to  Ga.As. 


Figure  1 shows  the  calculated  values  of  the  low-field 
(100  V cm)  drift  moliilitv  for  the  Ga|_.In,P,  .As,  qua- 
ternary alloys  with  a doping  level  of  lO’’'  cm'"  which 
are  lattice  matched  to  GaAs  and  InP.  These  mobilit\ 
values  are  obtained  by  a maximum  likelihood  estimation 
of  diffusion  coefficient,  and  by  using  the  Einstein  rela- 
tion lietween  diffusion  coefficient  and  mobilitv. ' ' This 
method  is  similar  to  a recent  technique  developed  bv 
Canali  c/  n/.'*"  to  calculate  low-field  drift  mobility,  and 
has  been  found  to  overcome  some  of  the  prolilems  in- 
lierent  in  the  application  of  the  Monte  Carlo  method  at 
low  electric  fields.  As  can  be  seen  in  Fig,  1,  the 
drift  moliility  of  the  alloy  matched  to  InP  lies  well  above 


TABI.K  I.  I'uinparl.son  of  lu.itci'iiil  iiiii'iiinctcrs  fm  :i  doping  iuvcl  o|  lid'  cm"  olitnincd  Irum  .Monte  t'arlo  ci'lculnlion.s. 


•wp  drift  mobility  .it  100  V cm. 


m.iximum  drift  velocity. 


' ImnijI  in.Tgnihide  of  iiKiximum  negative  dirierential  moliility. 
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the  value  for  GaAs  for  v ^ 0.  30  (a  < 0.  86).  Since  larRe 
values  of  electron  drift  mobility  are  desirable  for 
MESPET’s,  this  ransje  of  alloy  composition  looks 
especially  attractive. 

Velocity-field  calculations  have  been  made  alonn  both 
the  GaAs  and  InP  lattice-matched  lines.  For  the  GaAs 
substrate,  the  peak  velocity  for  the  quaternary  was 
found  to  continually  decrease  from  that  for  pure  GaAs. 
For  the  InP  substrate,  the  peak  velocity  was  found  to 
reach  a maximum  for  quaternary  alloy  with  v~0.6. 

Fijjure  2 shows  the  velocity-field  relationships  for 
the  quaternary  alloys  matched  to  InP  which  have  the 
larpest  calculated  values  of  peak  drift  velocity.  The 
value  for  a = 0.73,  y --0.6  of  2.95x10^  cm  'sec  is  the 
largest  we  have  calculated  for  any  material  with  a l)and 
gap  and  intervalley  separation  large  enough  to  l)c  useful 
for  MESFET's.  Also,  it  must  be  pointed  out  that  the 
values  presented  in  Fig.  2 assume  a completely  random 
alloy  with  no  sublattice  ordering,  so  that  the 

amount  of  alloy  scattering  in  the  Monte  Carlo  simulation 
is  maximum.  If  the  material  has  some  degree  of  order- 
ing so  that  the  amount  of  alloy  scattering  is  reduced, 
the  peak  velocities  could  be  larger  by  as  much  as 
10-20'T. 

Figure  3 illustrates  the  influence  of  alloy  scattering 
on  the  velocity-field  characteristics  of  the  quaternary 
with  v = 0.  73.  v = 0.6,  and  a doping  level  of  10'’  cm"’. 
Alloy  scattering  in  a random  alloy  can  be  seen  to  de- 
crease the  peak  velocity  and  increase  the  thre.shold 
field.  Also  shown  in  this  figure  are  the  velocity-field 
relationships  obtained  by  the  Monte  Carlo  method  for 
GaAs,  InP,  and  the  ternary  Ga„  jIn,,  jAs,  which  has 
been  proposed  as  a promising  material  for 
MESFET’s. Finally,  Table  I summarizes  the  mate- 
rial properties  of  interest  for  microwave  devices  for 
the  Ga„_j.,In„_73P,|  ,,AS|,_„  quaternary  alloy,  with  and  with- 
out random  potential  alloy  scattering,  along  with  the 
properties  for  GaAs,  InP,  and  Ga,,  .,In„  ^As.  The  ternary 
Ga„  5ln„  jAs  can  also  be  appro.ximatcly  lattice  matched 
to  InP  and  has  a band  gap  of  0.89  eV 

In  summary,  this  work  indicates  tliat  the 
Ga|.,In,P,..As7  quaternary  alloy  should  h.ive  desirable 
material  properties  for  microwave  devices.  The  cal- 
culated peak  velocity  for  the  quaternarv  system  lattice 
matched  to  InP  is  larger  than  the  calculated  value  for 


tlie  binary  or  ternary  III-V  materials.  Also  the  quater- 
naries have  a large  negative  mobility  and  a large  peak- 
to-valley  ratio,  which  should  be  useful  in  tran.sferred 
electron  devices.  The  work  points  out  a need  for  ex- 
perimental verification  of  the  results,  as  well  as  point- 
ing out  the  need  for  experimental  studies  of  ordering 
in  ternary  and  quaternarv  alloys,  since  the  actual 
amount  of  ordering  in  these  materials  will  determine 
which  will  be  most  useful  for  specific  applications. 
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Qualitative  Analysis  of  Thin  Gallium  Nitride  Films  with 
Secondary  Ion  Mass  Spectrometry 

J.  Edward  Andrews,*'  A.  P.  Duhamel,^  and  Michael  A.  Littlejohn 

Department  ot  Electrical  Engineering,  North  Carolina  Stale  University.  Raleigh.  North  Carolina  27607 


Gallium  nitride  (GaN)  thin  films  grown  from  the  vapor  phase 
by  pyrofyzing  GajCjHjjj'NHj,  were  analyzed  using  Secondary 
Ion  Mass  Spectroscopy  (SIMS).  Comparative  mass  spectra 
were  obtained  from  GaN  grown  by  two  other  laboratories  using 
two  different  techniques  and  were  found  to  be  similar.  The 
SIMS  technique  Identified  the  presence  of  Ga,  N,  and  O,  In 
all  the  GaN  films.  C was  delected  In  the  GaN  film  prepared 
by  pyrolysis.  Evidence  of  H as  part  of  the  Ionic  structure  was 
in  the  higher  mass  spectra  of  all  the  samples;  however.  It  corid 
not  be  determined  with  certalnity  If  the  source  of  H original  ed 
In  the  SIMS  Instrument,  the  samples,  or  both. 


Thin  gallium  nitride  (GaN)  films  grown  on  o-A1..0-i  sub- 
strates by  pyrolytic  decomposition  of  GafC^HjlrNH,-,  in  a 
chemical  vapor  deposition  system  have  been  descrilx'd  in  an 
earlier  publication  (1)  The  as-grown  films  were  typically 
yellowish  brown  in  color  instead  of  transparent  as  would  be 
eipected  for  a material  with  a .T  S-eV  energy  band  gap  (2l. 
This  yellowish  color  was  characteristic  of  the  films  grown 

* Present  address.  Hesearch  Triangle  Institute.  P O,  Hoi  121tM. 
Research  Triangle  Park.  N C 27709 

’Present  address,  I'S  Knergv  Research  and  Develi'imient  ^d 
ministration,  Washington.  DC  20'*-l.'» 


throughoitt  the  temperature  range  from  5(X)  to  1000  °C.  Since 
this  was  the  first  lime  that  this  growth  process  and  system 
had  been  ased  to  grow  GaN,  it  was  not  at  all  certain  that  the 
dejaisits  obtained  were  GaN,  especially  in  view  of  the  coloring 
noted  above.  Tbereftire  effort  during  the  early  part  of  the 
research  was  directed  toward  obtaining  a qualitative  analysis 
of  the  deposits  Several  anaKdical  techniques  were  eventually 
used  (electron  microprobe,  x-ray  diffraction,  electron  dif- 
fraction. secondary  ion  ma.ss  spectrometry)  with  each  pro- 
viding .some  useful  but  not  complete  information  concerning 
the  nature  of  the  deposits.  The  analytical  technique  reported 
here  is  worthy  of  note  liecaase  of  its  ability  to  delect  nilrrigen 
as  well  as  the  lighter  elements,  including  hydrogen,  and  was 
particularly  valuable  in  analyzing  deposits  that  were  tw 
amorphous  to  give  useful  diffraction  data. 

This  paper  reports  the  results  of  the  qualitative  analysis 
of  the  GaN  thin  films  using  the  technique  of  Secondary  Ion 
Mass  Spectrometry  (SIMS).  The  combination  of  sputter 
etching  and  mass  spectrometry  used  in  the  SIMS  technique 
permits  a convenient  and  rapid  qualitative  analysis  of  thin 
films  and  surfaces.  In  depth  discussions  of  the  SIMS 
technique  and  available  instrumentation  have  recently  been 
presented  (J  5).  A comparison  of  the  performance  of  the 
!->lMS  technique  with  Auger  Electron  Si>e,-lroacopy  (AES)  and 
X ra\  iiholiK-lectron  spectroscopy  (XPS)  given  in  Ref.  6 in- 
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Table  I.  Experimental  Conditions  Csed 
for  Sample  Analysii 

(1 ) Primary  ion  beam  gas:  argon” 

(2)  Argon  gas  pressure;  3x  10  * Torr 

(3)  Accelerating  voltage  for  argon  ions;  6 KeV 

(4)  Ion-beam  current  at  sample  surface;  1 50  >i  A 

” Oxygen  was  substituted  for  argon  and  mass  spectrum 
retaken  on  the  GaAs  sample  after  the  various  samples  had 
been  analyzed. 


10  15  20  25  JO  15  20  25 
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Figure  1.  Two  consecutive  mass  scans  of  sample  55S 

dicates  that  the  SIMS  technique  is  capable  of  a much  lower 
detection  threshold. 

For  purposes  of  comparison,  OaN  samples  furnished  by 
courtesy  of  IBM  and  RCA  were  al.so  an.alyzryl.  The  IBM  CaN 
was  prepared  by  sputtering  Oa  in  the  presence  of  N;  (7,  8). 
The  RCA  OaN  was  prepared  by  the  ammonolysis  of  GaCl  (9l. 

EXPKRIMKNTAL 

The  analysis  of  the  OaN  films  (plus  a sample  of  single  crystal 
undoped  Oa/Vs)  was  performed  using  an  li>n  Beam  Surfai  e Mass 
Analyzer  (ISMA).  This  instrument  and  it-s  operation  have  been 
described  in  detail  in  (.5).  The  analysis  of  the  undoped  GaAs  was 
performed  in  order  to  eliminate  possible  background  effects 
obtained  from  the  instrument  and  to  assist  in  the  interpretation 
of  the  spectral  data  obtained  from  the  OaN  films. 

The  analysis  of  the  OaN  film.s  and  purr  GaAs  were  performed 
under  constant  experimental  conditions  listed  in  Table  1 
Sample  chamber  pressures  prior  to  sample  analysis  were  or- 
idinarily  of  the  order  of  8 X 10  ''  Torr.  The  mass  resolution  of 
the  Ih.MA  was  typirally  0 I amir 

Mass  spectra  for  the  OaN  and  GaAs  were  obtained  over  the 
mass  ranges  10  25.  62-77,  and  80  95.  F.ach  range  was  scanned 
in  a 2-min  period  at  a rale  of  0.125  amu/s  with  several  sc.ans 
performed  con.secutivelv  t/>  obtain  data  on  the  variation  of  various 
peak  intensities  as  a function  of  time. 

Species  of  potential  interest  in  the  mass  ranges  10  25,  62  77, 
and  80  95  were  'T!.  “N.  'X);  *Gh,  '’On;  '’“Ga'^.  *Ga'‘N.  '*'Ga"'0. 
”Ga'’C,  ’'Oa"N,  and  ’'rta'*’0,  respectively. 

RESUl.TS  AND  IJISCl'SSION 
Figure  1 shows  spectra  resulting  from  the  first  two  mass 
.scans  of  OaN  sample  .5,5-.S  grown  at  6(X1  °C  by  the  process 
described  in  Ref.  1.  (The  ordinate  gives  a measure  of  the  ion 
br-am  current,  in  A,  prixliiced  by  the  serundarv  ions  sputtererl 
from  the  surface  of  the  sample.)  Figure  1 indicates  hydr's 
carbon  fragments  tieing  detected  only  in  the  first  few  minutes 
of  analysis.  Normally,  hydrocarlam  fragment.s  detected  with 
this  technique  and  exhibiting  this  Fx'hnvior  (disapix-aring  after 
a few  minutes  of  etrhingl  are  usually  attributed  to  vacuum 


MASS /CHARGE  (omu) 
Figure  2.  Miiss  spectrum  ol  GaAs 


system  pump  oil  that  has  accumulated  on  the  surface. 
However,  0a.A.s  spei  lra  obtained  under  identical  experimental 
conditions  (.see  F'ijttire  2)  showed  no  indication  of  hydrexarbon 
fragments.  This  would  strongly  suggest  that  the  hydrixarbon 
fragments  observed  in  the  OaN  spectra  were  not  an  in.stru- 
ment  artifact. 


The  sample  55-S  mass  peaks  observed  at  14  and  16  amu 
[x-rsisted  with  continued  etching  while  those  at  12  and  13  amu 
gradually  disappeared.  The  peaks  at  12,  13,  and  16  amu  were 
interpreted  as  ('■Cl",  (CH)",  and  respectively. 

The  interpretation  of  the  peak  at  14  amu,  however,  is  more 
complicated  since  it  could  be  due  to  (CHo]",  I'^N]  or  both. 
The  gradual  disappearance  of  the  [’■C)"  and  [CH]"  peaks 
suggests  that  the  related  [CH>]"  peak  should  similarly  dis- 
apixar.  Since  the  ix’ak  at  14  amu  [x*rsLsted,  it  could  eventually 
be  interpreted  as  being  primarily  ('*N]"  ns  indicated  in  the 
second  scan  in  Figure  1.  'Fhe  peak  at  14  amu  was  seen  to  still 
be  very  evident  in  F'igure  3 which  is  the  spectrum  after  8 to 
10  min  of  etching.  The  absence  of  similar  peaks  in  the  GaAs 
sjxx-tra  (Figure  2)  indicaU-s  that  background  interferences  were 
not  responsible  for  the  observed  peaks  in  the  OaN  spectra. 

It  was  noted  (compare  F'igures  2 and  3)  incidentally  that 
the  On  ion  yield  for  sample  55S  was  more  than  1 order  of 
magnitude  larger  than  that  of  the  GaAs  sample.  Sample  5.5S 
was  determined  previously  to  be  nearly  amorjihous  in  nature 
and  might  therefore  he  expected  to  yield  Oa  ions  much  more 
easily  than  the  single  crystal  GaAs  used  for  the  reference. 

F'igure  3 shows  the  mass  spectra  obtained  for  sample  .55S 
through  the  ma.ss  range  80  95  amu.  No  significant  mass  peaks 
were  observed  for  the  GaAs  in  this  mass  range.  The  55S 
spectra,  however,  show  distinct  peaks  at  85,  86,  87,  88,  and 
89  amu.  Weak  peaks  were  also  observed  at  83  and  84. 

It  is  interesting  here  to  compare  these  data  with  the  GaAs 
ma.ss  sp<>ctrum  obtained  in  the  80  95  amu  range  asing  an  O;* 
primary  ion  beam  (see  F'igure  4).  As  was  seen  in  F'igure  2, 
the  Oa.As  spectrum  (sputtered  with  Ar)  gave  no  observable 
peaks  in  the  80  90  amu  range  whereas  F'igure  4 shows  peaks 
that  were  similar  to  the  OaN  sixetra  (the  oxygen-containing 
fragment-s  were  formed  as  a result  of  ion-molecule  reactions 
between  the  sample  and  O;"  ions  from  the  beam).  This 
comparison  would  strongly  indicate  that  the  OaN  films 
contained  oxygen 

F'lgurf  4 also  indicates  the  existence  of  hydrogen  in  some 
of  the  fragments  detected  Roth  (AsO]*  and  lAsOH)' 
fragment.s  were  Ixlieved  to  lx  responsible  for  the  peaks 
observed  at  91  and  92  amu.  respectively  The  Oa  from  the 
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Figure  3.  Mass  spectrum  of  sample  55S  (GaN)  with  the  mass  range  80-95  amu  shown  at  two  dittcjfent  sensitivites  (this  spectrum  persisted  throughout 
the  remainder  of  the  analysis  of  this  work) 
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Figure  4.  Mass  spectrum  of  GaAs  sputtered  by  O;*  primary  ion  beam 
Evidence  of  formation  of  both  gallium  and  arsenic  oxide  fragments  in 
the  analytical  environment 


It  cannot  bo  stated  at  this  timr  whothrr  the  hydrogen  is 
rontrdtuted  hv  the  IS.MA  environment  or  whether  it  is  an 
impurity  common  to  both  the  (iaAs  and  GaN  samples.  A 
fragment  corresponding  to  (GaCI*  IGaCHJ*  was  not  de- 
tected at  81  and  82  amu.  respectively,  and  therefore  cannot 
make  any  delectable  contribution  to  the  83.  84  amu  peaks. 

The  likelihood  of  |‘^'‘GaGH..|‘  and  |'  'GaCH  ,)*  was  con- 
sidered in  the  interpretation  of  the  (leaks  considered  at  8.3  and 
84  amu  since  it  is  here  that  |' 'GaN|*,  and  I'^^'GaNHj*  wmild 
be  detected  if  they  exist  and  could  result  in  an  ambiguity  in 
interpretation  The  identificBlion  of  the  fragment  ions  at  8.3 
and  84  amu  as  |*''GaNl*  and  |''"f;aNHl*.  respectively,  is 
supported  by  comparing  the  mass  spectia  of  (iaN  samples 
grown  by  the  two  other  pnaesses  referenced  earlier. 

Comparison  of  Mass  Spectra  for  Ga.V  Grown  by  Other 
Techniques.  Two  other  samples  of  GaN  grown  by  the  two 
techniques  referenced  earlier  (idenlifii'd  as  the  IBM  and  RCA 
.samples)  were  analyzed  while  still  using  the  Ga.As  for  a 
reference  The  sjiectra  for  the  IBM  sample  are  shown  in 


Figure  5 while  the  spectra  for  the  RCA  sample  are  shown  in 
Figure  6. 

These  GaN  .samples  resulted  in  spectra  similar  to  5bS 
throughout  the  83-89  amu  range.  In  addition,  these  spectra 
were  also  simil.ar  to  the  one  shown  in  Figure  4 (Ga.As  sputtered 
with  oxygen)  in  the  85  89  amu  range.  The  RCA  sample  was 
small  and  thus  exposed  the  edge  of  the  sapphire  (o-.MO.il 
sulistrate  to  the  primary  ion  lieam;  thas  the  oxygen  indicated 
by  the  spectra  for  the  RCA  sample  (i.e..  through  (GaO)* 
fragments)  could  easily  have  resulted  frvim  the  ion  molecule 
reaction  between  sulistrate  furnished  oxygen  and  sample 
supplied  Ga  atoms.  The  RCA  sample  contributed  rather 
definite  peaks  at  83  and  84  amu  which  were  interpreted  at 
('■‘GaN)*  and  |‘'‘‘GaNHl',  respectively,  which  supports  a 
similar  interpretation  in  the  55S  spectra. 

.Sample  55S  and  the  IBM  sample  did  not  experience  any 
substrate  interference;  therefore  oxygen  indicated  in  their 
spectra  was  interpreted  as  originating  from  the  s.smpies. 

As  stated  earlier,  hydrogen  is  also  very  much  in  evidence 
in  these  spectra,  hut  cannot  be  identified  as  originating  from 
the  sample  or  the  ISMA. 

CONCI.liSlON.S 

The  qualitative  analysis  of  deposits  obtained  from  the 
[lyrolysis  of  Ga(C  ,H  ,)  .-NH;,  has  shown  that  the  sample 
contains  gallium,  nitrogen,  oxygen,  and  c.arlxin,  the  latter  licing 
concentrated  near  the  surface  of  the  sample.  While  x-ray 
diffraction  was  relied  on  for  the  verification  of  the  deposits 
to  he  GaN,  it  is  significant  that  SIMS  provided  a means  of 
directly  detecting  the  presence  of  nitrogen  in  the  deposits 
HydriKarhons  from  the  vaniiim  system  oil  were  ruled  out  as 
a source  of  carlion  because  of  the  absence  of  hvdrocarhons 
in  the  Ga.\s  sfiecira  HydriK-arbi'iis  detectr-d  in  the  S-5,5S  film 
would  la’  significant  since  Ihi.x  was  suspected  to  lie  the  rea.son 
for  the  yellow- brown  rolor  that  has  Ix’en  characteristic  of  the 
films  grown  from  Gn(t’ .H-,),i-NH.i.  It  was  interesting  to  note 
that  samjile  ,5,'iS  has  changed  color  (in  the  region  lieing  an- 
alyzed) from  vellowish  brown  to  gray  No  explanation  ran 
he  offered  for  this  rhange  in  rolor  at  the  present  time 

The  GaAs  sample  permitted  establishing  a background 
baseline  against  which  the  GaN  spectra  could  be  compared. 
The  ma.s.s  sjiectra  obtained  from  the  (iaN  grown  by  two  olhet 
lalviralories  adder!  ixinfidenre  to  the  interjiretation  of  the  nia-s.*- 
pr-aks  as  did  the  s|K’<tra  obtained  by  using  the  oxygen  primary 
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Figure  5.  Mass  spectrum  of  IBM  C^N 
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Figure  6.  Mass  spectrum  of  RCA  GaN 

ion  beam  on  the  OaAs  sample  at  the  end  of  the  analysis. 

Some  of  the  fragment  ions  occurred  at  an  m/c  that  could 
be  explained  tmly  through  the  addition  of  one  or  two  hydrogen 
atoms  (or  ions).  Although  the  analytical  environment  of  the 
ISMA  could  itself  be  a source  of  hydrogen,  one  cannot  rule 
out  the  actual  growth  process  and  systems  as  a source,  es- 
pecially where  reactants  such  as  NHi  or  Cia(C.H,s)i-NH  , arc 
involved.  It  was  in  fact  the  hydrogen  question  (hat  motivated 
the  analysis  of  the  IBM  sample  which  did  not  intentionally 
use  any  hydrogen-containing  reactants.  Only  extremelv  careful 
experimental  procedure  ran  hope  to  pinpoint  the  sourrels) 
of  the  hydrogen.  It  is  emphasized  here  that  the  evidence  for 
hydrogen  persisted  even  after  extended  lengths  of  s|)uttrr 
etching 


I'pon  reviewing  these  spectra,  it  appears  unnecessary  to 
seek  (GaN)*  ions  directly  as  evidence  of  the  N in  GaN  for  the 
experimental  conditions  desrrilied  here.  Nitrogen  ions  de- 
tected at  14  amu  are  attributed  to  the  sample  and  not  to  the 
instrument  background  for  the  conditions  of  our  analysis. 
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Velocity-field  characteristics  of  GaAs  with 
conduction-band  ordering^^ 
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This  paper  describes  Monte  Carlo  calculations  of  vclocity-fteld  characteristics  for  GaAs  using  the  rccenl 
esperiinental  conduction-band  ordering  of  Aspnes,  which  places  the  /.((III)  conduction-band  minima 
lower  in  energy  than  the  TJ(IOO)  minima.  These  calculations  use  intervalley  deformation  potentials  which 
give  the  best  fit  to  recent  high-ficId  drift  velocity  measurcnienls,  and  at  the  same  time  give  gOtxJ 
agreement  with  accepted  peaV  velocity  and  threshold  field  values. 

PACS  numbers  71  25.Tn,  72  20  Ht 


The  recent  Schottky -barrier  electroreflectance  mea- 
surements using  synchrotron  radiation'’^  have  offered 
evidence  that  the  Ljllll)  conduction  -band  minima  in 
GaAs  is  lower  in  energy  than  the  ,V-{IOO)  minima.  Tliis 
is  contrary  to  previous  conduction -band  ordering,^  and 
the  observation  has  led  to  a great  deal  of  activity  toward 
evaluating  consequences  of  it.  Hydrostatic  pressure 
and  uniaxial  stress  measurements  in  GaAs  *•*  suggest 
the  rf-Lj-A'I  ordering  and  are  in  agreement  with  Monte 
Carlo  calculations  based  on  this  assumed  band  struc- 
ture. Separate  hydrostatic  pressure  measurements' 
and  uniaxial  stress  measurements^  complement  the 
above  experiments,  and  the  results  of  resonant  Kaman 


scattering  have  also  been  interpreted  in  terms  of  this 
conduction -band  structure.*  It  appears  that  the  evidence 
being  accumulated  indicates  that  the  rf-Z.|-A'|  conduc- 
tion-band ordering  is  correct  for  GaAs. 

There  are  Monte  Carlo  calculations  of  velocity -field 
characteristics  for  GaAs  using  the  previously 
accepted  conduction -band  ordering  which  are  in  good 
agreement  with  experimental  data.  These  calculations 
are  based  not  only  on  the  band  structure,  but  require 
a knowledge  of  material  parameters  which  have  not 
been  measured  and  which  are  used  as  adjustable  param- 
eters in  obtaining  agreement  between  these  computer 


TABLE  I.  Ga.As  material  p.aranieters  ii.-ied  in  the  calculations. 


Bulk  material  paranu’ters 

Parametor  Value 

Parameter 

Value 

L.itticc  con.st.'uit  5.  fi-12  x ur* 

IX)  phonon  enorf^y 

0.  03536  eV 

Density  S.Hfigcm’ 

Sound  Vflocit\’ 

5.  24  X]  0*  cm/  sec 

F lectron  affinity*  4.07  eV 

Optical  dielectric  constant 

10.  92 

Plc7oc*1ectric  constant  0.16  C 

Static  dielectric  constant 

12.90 

V'alley-depcndent  material  parameters 

Parameter 

r (000) 

Conductum-band  valley 

X(IOO) 

/.(111) 

Acoustic  deformation  potential  (eV) 

7.0 

9.27 

9.2 

Effective  mass  (rn'/nif) 

0.003 

0.58 

0.222 

Nonpnralx)l(city  (eV'"') 

0,610 

0.20t 

0.461 

Energy  b;md  (eV) 

l.JSS 

1 . 001 

1.  709 

(relative  to  valence  band) 

Optical  deformation  itotcntial  (eV/cm) 

0 

0 

.3  X \ o" 

Optical  phonon  energy  (cVj 

0.n,'t43 

Intervalley  deformation  jKiU  ntial  (rV/cm) 

from  r(OOO) 

0 

1 Xio' 

1 xio’’ 

from  xaoo) 

1 xio' 

7X]0' 

5X10* 

frt>ni  LO-ll) 

1 Xio’ 

5 xio' 

1 xio’ 

Intervalley  phonon  ener^  (eV) 

from  noOO) 

0 

0.n29<J 

0.0278 

from  XOOO) 

0.0259 

0.0299 

0.0293 

from  L(X\\} 

0.0278 

0.0293 

0.  0290 

Numl.)er  of  equivalent  valleys 

1 

3 

4 

•’This  work  was  sunwrted  by  a research  grant  from  tlie  (iffice 
of  Naval  Research,  Washington,  D.C. 
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FIG.  1.  Drift  velocitj'  versus  electric  field  intensity  for  Ga-^s 
with  a doping  level  of  lO'®  cm"®.  The  solid  curves  are  for 
T-L-X  conduction-band  ordering  with  c\'  and 

= 0.192  eV.  The  dotted  curve  is  for  T -X-L  ordering 
using  Fj- = 1.439  eV,  £;f  = l.bG9eV,  £i  = 2.0eV,  andwi'Ot) 

= 0. 39m(|.  The  b.ars  show  the  range  of  repoi  ted  experimental 
values.  The  values  of  are  2.5x10*  eV/cm  (curve  A), 

5X10*  eV/cm  (curve  B),  1 xio®  eV/cm  (curve  C),  2.5 
xi0*eV/cni  (curve  D). 

simulations  and  the  experimental  velocity -field  curves. 
The  question  arises  as  to  what  effect  the  new  conduc- 
tion-band ordering  has  on  the  Monte  Carlo  calculations, 
and  what  material  parameter  adjustments  have  to  be 
made  in  order  to  obtain  the  same  good  agreement  with 
experiment.  Previous  authors*'®  have  addressed  this 
question  to  some  extent,  although  they  have  been  pri- 
marily concerned  with  pressure  and  stress  effects  and 
with  electric  field  intensities  near  the  threshold  field, 
i.e. , £s5  kV.'cm.  In  addition,  much  experimental 
velocity -field  data  taken  to  dale  has  been  restricted  to 
electric  field  values  below  about  20  kV  cm.  Data 
taken  at  higher  electric  fields  (20  — 100  kV-  cm)  show- 
much  variation  in  absolute  values  between  different 
authors,  probably  due  to  differing  experimental 
tecliniques'®  and  possibly  due  to  contact  effects.®  This 
makes  theoretical  calculations  in  the  high-field  region 
uncertain. 

Houston  and  Evans'®  have  very  recently  used  the  tinie- 
of -flight  technique  to  obtain  what  appears  to  be  the  most 
reliable  velocity -field  data  taken  to  date  for  fields  be- 
tween 20  and  100  kV  cm.  Using  this  data,  and  other 
reliable  low -field  data  near  the  threshold  field,  as 
an  experimental  basis  this  paper  describes  Monte  Carlo 
calculations  for  GaAs  from  low  fields  to  100  kV  'cm 
considering  the  T\-Ll-X\  conduction -band  ordering. 
Intervalley  deformation  potentials  have  been  used  which 
give  reasonable  fits  to  high-field  data  and  which  also 
give  good  agreement  with  data  near  the  threshold  field. 

The  Monte  Carlo  transport  program  has  beoii  pre- 
viously described.'*''*  and  is  very  similar  to  others 
reported  for  CaAs  studies.*""'"  Table  1 summarires 
the  material  parameters  used  to  calculate  velocity -field 
characteristics  for  GaAs  with  T-T-.V  conduction -b.and 
ordering  in  our  work.  These  (laramoters  are  essentially 
the  same  as  those  employed  in  other  published  siimila- 
tions,’""'"  except  for  the  L and  X band  gaps  and  a few 
other  changes  to  be  discussed  in  the  text. 


Figure  1 shows  velocity -field  curves  for  GaAs  calcu- 
lated with  T-L-X  conduction-band  ordering.  The  energy 
band  separations  of  0.  284  eV  and  i.  0.  192 

eV  originally  reported  by  Aspnes,'"  llie  ,Y  effective  mass 
of  0.  39<N|),"  and  the  equivalent  intervalley  deformation 
potential  in  the  L valley,  Dt-o  3x10*  eV/cm  " have 
been  used.  Sliown  for  comparison  is  a curve  computed 
using  r-A'-£  ordering  with  the  previously  accepted 
material  parameters."  The  T-L  intervalley  deforma- 
tion potential,  has  been  varied  from  2.5x10*  to 

2.5x10*  eV,-cni.  Besides  the  variations  in  tlie 

effective  mass  )n*(.'.),  the  deformation  potential 
and  the  different  conduction -band  ordering  and  separa- 
tions, all  other  m.aterial  parameters  have  the  value.s 
given  in  Table  1.  Figure  1 also  shows  the  range  of  ex- 
perimental values  of  drift  velocity  for  fields  between 
2 and  10  kV  cin."''®  Ionized  impurity  scattering  has 
been  included  in  the  simulations  with  an  impurity  con- 
centration of  lO'*  cm'®.  The  iiiclu.sion  of  this  scattering 
mechanism  is  kniv  n to  reduce  the  ma.ximum  velocity 
and  to  increase'  me  threshold  field'"-'®''*  and  was  not 
included  in  previous  pressure  effect  studies.  *'*  It  is  to 
be  noted  that  for  these  energy  separations  the  .agher 
values  of  (~  10*  eV  cm)  are  required  to  move  the 
calculated  curves  into  the  range  of  experimental  values 
near  the  th-eshold  field.  The  low-field  (100  \'  cm) 
drift  mobility  for  the  curves  of  Fig.  1 has  been  computed 
by  the  maximum -likelihood  estimation  technique''*''*  and 
is  6600  cm®  V sec. 

Our  calculatioi.s  indicate  that  the  only  material 
parameter  which  can  be  reasonably  varied  to  extend  the 
velocity-field  curves  toward  those  for  the  T-X-L  order- 
ing is  the  energy  separation  AFp.j^.  Other  authors*'® 
have  had  to  use  F-i.  .separations  larger  than  0.  284  oV 
to  obtain  agreement  between  stress  experiments  and 
Monte  Carlo  calculations.  Very  recently,  .Aspnes  has 
used  uniaxial  stress  experiments  to  infer  a F-L  separa- 
tion of  A/.'r.i  - 0.  330  1 0.  05  eV.  ‘ Our  calculations  have 
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been  repeated  for  higher  separations  ajid  Kig.  2 shows 
the  results  (or  AEp-i  = 0- 33  eV  and  eV. 

The  velocity  near  the  threshold  field  is  very  similar  to 
both  experimental  data  and  previous  simulations.  How- 
ever the  slope  of  the  velocity -field  curve  above  the 
threshold  field  appears  to  be  too  small.  For  AEy--l 
= 0.38  eV  and  = 0.192  eV  the  peak  velocities 

are  somewhat  higher  (t’„u  = 2. 1 x lO’  cm/sec  for  £=4.7 
kV/cm  with  = I x 10*  eV/cm)  and  the  slope  above 
threshold  is  smaller  still.  All  our  calculations  indicate 
that  A£r.i,  » 0.  33  eV  and  Or-t  ^ 1 ^ 10*  eV/cm  are  re- 
quired to  give  agreement  between  the  Monte  Carlo 
simulations  and  experimental  velocity -field  data  near  the 
threshold  field. 

The  value  of  the  equivalent  LI  intervalley  deforma- 
tion potential,  Oi-r,  w-as  varied  to  change  the  slope  of 
the  velocity -field  curve  above  threshold.  It  was  found 
that  as  Dfi  varied  from  3xio*  to  above  10*  eV  cm, 
the  velocity  at  16  kV/cm  could  be  decreased  while  the 
velocity  at  5 kV'  cm  was  effectively  unchanged,  indicat- 
ing the  desired  increase  in  slope.  The  be.st  value  of 
^L-L  seems  to  be  about  (1— 2)xi0*  eV/cm  which  results 
in  a value  of  velocity  at  16  kV  cm  equal  to  the  value 
required  by  the  experimental  data  and  the  previous 
T-X-L  curves  shown  in  Figs.  1 and  2.  The  actual  best 
value  of  Di.i  depends  somewhat  on  the  AFp.j,  separation, 
and  good  fits  to  the  data  can  be  achieved  for  fields  up 
to  20  kV  cm  using  = 10*  eV  cm  and  33 

eV. 

The  value  of  - 10*  eV  cm  also  gives  reasonable 
agreement  with  the  best  high -field  data  available.  Fig- 
ure 3 shows  this  data  compared  with  the  Monte  Carlo 
calculations  using  aF^.^,  = 0.  33  and  0.  38  eV.  Two  sets 
of  data  are  shown  representing  original  measurements'* 
and  a very  recent  refinement'*  both  taken  by  the  time- 
of-flight  technique.  The  Monte  Carlo  data  for  AFp.j, 

- 0.  33  eV  results  in  a nearly  correct  field  dependence 
for  the  velocity.  However,  because  of  the  higher  calcu- 
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FIG.  a.  Drift  velocity  versus  electric  field  intensity  for  Gri-As 
with  a doping  level  of  10*'  ctn**.  The  dotted  curves  arc  exiieri- 
nicnlal  data  (Curve  1 — le.ist-squ.nres  fit  to  d.iLi  of  Kef.  19 
curve  2— <1ati  of  Het.  13).  The  solid  curves  are  Monte  Carlo 
data  (curve  A — Afr.^  0.33  cV,  iu*iA')  0.3''ii,j  curve  li — 
AFj-.,  = 0.  33  eV,  » • (Jf)  O.Ssni,  lurveC — A/.,-.j^  0.3s  cV, 

m*  (A3  » 0.  ."iRm,,  For  .ill  curves  AF^.j  0.102  eV). 
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FIG.  4.  Drift  velocity  versus  temperature  for  GaAs  with  a 
doping  level  of  10*'  cm"’’,  E 100  kV  'em,  AJTj..^  = 0. 33  eV,  ."uid 
A£^_jf  --  0. 1S2  eV.  Curve  A — w’  tV)  0. 30m(|,  curve  B — m'  tY) 

- O.TSnij,  cuive  C — m‘  tX')  - 0.  OSm^,  curve  D — -AHj— i ‘ 0-  33  eV, 
Af:^.^  = 0.142  eV,  m*(A3  = 0.39TOo. 


latad  values  compared  with  the  experimental  values, 
tiie  effective  mass  in  the  valley  was  increased  to 
reduce  the  high-field  velocity.  '*  The  results  of  this 
calculation  are  shown  in  Fig.  3.  The  /.|  mass  could  also 
be  varied  to  reduce  the  velocity,  taut  the  value  tor  this 
band  seems  mure  certain  than  the  mass  for  the  .X''f 
band,*  and  varying  the  /.?'  mass  would  affect  the  near- 
threshold drift  velocity.  A value  for  the  effective  mass 
ratio  in  the  .Y  valley  of  0.  58  gives  the  best  fit  to  the 
data  for  A£r-t  = 0- 33  eV,  while  increasing  the  separa- 
tion to  0.  38  eV'  causes  a deviation  at  the  lower  fields. 
This  agreement  is  also  manifested  in  the  temperature 
dependence  of  the  near -saturation  velocity.  Figure  4 
shows  the  measured  temperature  dependence  of  the 
drift  velocity  at  100  kV  cm  compared  to  the  Monte  Carlo 
calculations  as  the  mass  in  the  A’|  valley  is  varied.  The 
temperature  dependence  and  field  dependence  of  the 
high -field  drift  velocity  obtained  from  Monte  Carlo 
calculations  are  in  agreement  with  experiment  using 
m*(X'|)--  0.  SSm,. 

Another  signific.ant  \ . Triable  is  the  Lf-.V-  cuiergy 
spacing.  aE^^.,  . A.spncs  has  reported  this  parameter 
to  be  0.  170±0.  03  eV  at  2 K,*  with  a value  of  0.  192 
> 0.04  oV  '*  at  300  K.  We  have  studied  variations  of 
aE'^.jj  on  the  Monte  Carlo  calculations  for  GaAs  with 
r-E -,Y  conduction -band  ordering.  For  example,  curve 
D of  Fig.  4 shows  the  100-kV  cm  drift  velocity  versus 
temperature  for  A£r./,  = 0.'33  cY.  m*(.Y)-  0 39i(;(i,  and 
a£;_.j  0.  142  eV.  There  is  essentially  no  change  with 

Die  decrease  in  AEj^.i.  In  fact.  AE'^.jt  would  have  to  be 
decrea.sod  below  0. 1 eV  to  result  in  significant  changes 
ui  any  results  presented  so  far,  and  this  is  viewed  as 
an  unreasonable  value  fn  light  of  Aspnes’s  data.  At  the 
same  time,  other  authors*’*  have  used  AF^.,  as  small 
a.s  0.02  eV  in  interpreting  pressure  effect  data.  How- 
ever, there  is  a question  regarding  this  data.'  and 
the  present  author.s  view  such  small  values  of  A/.^.,  as 
somewhat  unrcasonaole  until  further  results  are 
achii’vod. 
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In  summary,  Monte  Carlo  calculations  and  experi- 
mental velocity -field  data  are  in  agreement  using  the 
r{-I,|-X|  conduction -band  ordering.  The  major  param- 
eter changes  required  from  those  used  in  previous 
T-X-L  band-ordering  calculations  to  achieve  this  agree- 
ment are  1 xlO’  eV/cm,  = 1 x 10*  eV/cm, 

and  m*(Xf)  = 0.  58mo,  with  a minimum  T~L  energy 
separation  of  0.  33  eV.  These  values  of  deformation 
potentials  do  not  coincide  with  the  values  = 2.8x10* 
eV/cm  and  0^.1,  = 1.  8 x 10*  eV/cm  required  to  give  the 
best  agreement  to  pressure  and  uniaxial  stress  experi- 
ments.* However,  the  actual  variation  and  amount  of 
increase  of  peak  velocity  and  threshold  field  for  hydro- 
static pressures  up  to  10  kbar  appears  to  be  in  ques- 
tion,* and  the  previous  Monte  Carlo  calculations  probably 
need  to  be  reevaluated  in  light  of  new  data.  * In  view  of 
this  uncertainty,  the  calculations  made  here  are 
presented  with  material  parameters  which  give  agree- 
ment with  experimental  data  for  electric  fields  as  large 
as  100  kV/cm.  Additional  calculations  as  well  as  ex- 
perimental verification  of  deformation  potential  values 
should  be  made. 
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Energy  band  gap  and  lattice  constant  contours  are  presented 
for  the  nine  quaternary  alloys  formed  from  Al,  Ga,  In  and 
P,  As,  Sb . The  quaternary  bandgaps  wei’e  obtained  using  an 
interpolation  formula  proposed  by  Moon  ^ al . The  quater- 
nary lattice  constants  were  obtained  by  use  of  a linear 
interpolation  technique  using  the  binary  lattice  constants 
as  boundary  values. 

Key  words:  quaternary  alloys,  bandgap,  lattice  constant. 

Introduction 

There  has  been  considerable  interest  in  the  quater- 
nary  III-V  semiconductor  materials  for  many  applications, 
such  as  electro-optics  and  microwave  devices  [l-6j.  One 
significant  reason  for  this  interest  is  the  ability  to 
synthesize  materials  with  a fixed  lattice  constant  and  a 
variable  range  of  energy  bandgaps  or,  conversely,  to 
synthesize  materials  with  a constant  energy  bandgap  and 
a variable  lattice  constant  [7]. 

*This  work  was  supported  by  a research  grant  from  the 
Office  of  Naval  Research,  Washington,  D.C. 
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In  order  to  predict  the  energy  baridgap  or  lattice 
constant  for  a quaternary  material  a two-step  sequence  is 
involved.  Initially,  it  is  necessary  to  have  available 
experimental  (or  theoretical)  values  of  these  material 
parameters  for  the  four  possible  binary  III-V  constituents 
of  the  quaternary  compound.  These  values  are  readily 
available  for  most  III-V  materials  in  the  literature  [8], 
There  are  also  four  possible  ternary  (pseudobinary)  III-V 
systems  in  a given  quaternary  compound.  In  the  first  step, 
the  energy  bandgap  and  lattice  constant  as  a function  of 
alloy  composition  for  each  ternary  system  is  computed 
from  the  known  binary  material  parameters.  This  computa- 
tion is  on  a good  experimental  and  theoretical  basis 
[9-13]  involving  the  use  of  Vegard's  law  and  the  concept 
of  the  "bowing  pai'ameter"  for  the  energy  bandgap  [14]. 

The  final  step  involves  the  use  of  the  compositional 
dependence  of  the  energy  bandgap  and  lattice  constant  for 
the  four  tern^iry  combinations  to  interpolate  to  the  desired 
quaternary  energy  bandgap  and  lattice  constant.  This  final 
interpolation  step  is  somewhat  empirical  in  nature.  There 
appear  to  be  several  such  interpolation  techniques  used  in 
the  literature  [6,7,15]  and  while  these  are  similar,  they 
also  have  non-trivial  differences.  This  paper  will  discuss 
this  critical  interpolation  procedure  in  view  of  the  small 
amount  of  experimental  quaternary  data,  and  will  present 
calculations  fur  energy  bandgaps  and  lattice  constants 
for  the  nine  quaternary  alloys  formed  from  AL,  Ga,  In  and 
P,  As,  Sb  using  one  of  these  interpolation  techniques. 


Notation 


There  appears  to  be  no  generally-accepted  scheme  for 
symbolically  describing  the  III-V  quaternary  alloys.  In 
this',  paper  the  following  convention  has  been  adopted. 
Within  both  the  group  III  and  group  V pairing,  the  first 
position  is  occupied  by  the  element  with  lowest  atomic 
numl.er.  The  composition  variables  are  associated  with  the 
four  elements  in  the  order  1-x,  x,  1-y,  y.  For  example,  a 
.quaternary  alloy  will  be  denoted  as  A,  B C,  D . Here  A 
and  B are  group  III  elements  with 
number  than  B.  Likewise,  C and  D 
C having  lower  atomic  number  than 
authors  [6,7],  a quaternary  alloy 
or  lattice  constant)  is  described 


A,  B C,  D . 

1-x  X l-y  v 
A liaving  lower  atomic 

are  group  V elements  with 
D.  Following  previous 
pai-ameter  (e.g.,  bandgap 
by  a surface  Q(x,y)  over 
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the  x,y  composition  plane  (0<x<l,  0_<y<l).  At  the 
corners  (x  and  y equal  to  zero  and/or  one)  the  values  of 
the  parameter  for  the  four  binary  elements  Q(0,0)  = B , 
Q(1,0)  = B , Q(l,l)  = Bg  and  Q(0,1)  = are  obtained. 
Along  the  boundaries  of  the  plane  the  parameter  for  the 
four  ternary  elements  Q(x,0)  = T (x),  Q(l,y)  = T (y), 
Q(x,l)  = T^gCx)  and  Q(0,y)  = obtained. 


Estijraj^ion  of  Alloy  Parameters 

Many  ternary  alloy  parameters  (e.g.  lattice  constant) 
are  obtained  by  linear  interpolation  from  those  of  the 
constituent  binary  compounds,  i.e., 

T.  .(x)  = xB.  + (l-x)B.  (1) 

For  some  parameters  the  theoretical  variation  with  compo- 
sition is  nonlinear  [14,16].  For  the  direct  bandgap  for 
example,  Thompson  and  Woolley  [9]  have  shown  that 


T.,(x)  = xB.  + (l-x)B.  - C..x(l-x)  (2) 

3 11] 

where  C.  . is  the  bowing  parameter  for  the  ternary  alloy 
bandgap  '^qj* 

In  the  absence  of  definitive  theories  for  quaternary 
parameters,  estimates  of  a quatei'nary  alloy  parameter 
Q(x,y)  must  be  obtained  by  interpolation  from  the  four 
ternary  alloy  parameters  T...  Various  interpolation 
schemes  have  been  proposed^'^  Onton  and  Chicotka  [15] 
used  the  solution  of  Laplace's  equation  subject  to  the 
boundary  conditions  Q(x,0)  = T^^(x),  etc.  In  Monte  Carlo 
transport  studies  [6]  the  present  authors  have  used  the 
interpolation  equation 


Q(x,y)= 


x(l-x)[  ( l-y)!^^^^  ) ]+y(  1-y  )[  ( l-x)T^^^(y  ) + xT22(y  ) ] 

x(l-x)+y( 1-y ) 


(3) 


This  interpolation  equation  reduces  to  the  ternary  param- 
eters on  the  quaternary  plane  boun/.aries  and  to  the  average 


3 


Glisson 


73 


of  the  ternary  parameters  at  the  midpoint  (x=0 . 5 ,y=0 , 5 ) 
of  the  compositional  plane.  This  is  incidentally  the 
solution  to  Laplace's  equation  if  the  T,.  are  linear  as 
in  Eq.  (1).  Moon  et  al.  [7]  have  proposed  similar  schemes 
for  the  estimation  of  lattice  constant  and  bandgap.  For 
the  lattice  constant,  Moon  gives  (in  the  above  notation) 

Q(x,y)  = B^+(B2-B^)x  + (B^-B^)y  + (B^-B2+B2-B^)xy  (4) 

For  the  lattice  constant,  the  T..  are  assumed  to  be 
linear  in  their  arguments,  so  the  soitition  to  Laplace's 
equation  is  identical  to  Eq.  (3),  which  in  turn  also 
reduces  to  Eq.  (4). 

For  the  bandgap.  Moon  et  have  used  the  equation 

Q(x,y)  = (l-x)T^^^(y)  + x T^^Cy)  - A (5) 

where  the  T.  . are  determined  from  Eq.  (1)  and 
11 

A = x(l-x)[(l-y)C^2+yC^^3]  + y(  1-y ) [ ( l-x)C^^+xC23  ] (6) 

An  important  difference  between  this  bandgap  estimate  and 
that  obtained  from  Eq.  (3)  is  in  the  manner  in  which  the 
ternai’y  bowing  parameters  C..  enter  the  calculation.  In 
Eq.  (3),  the  ternary  bowing^^tfects  are  included  in  the 
T. whereas  in  Moon's  method  the  bowing  enters  as  a 
separate  quaternary  bowing  parameter  term  given  in  Eq.  (6). 
At  the  center  of  the  composition  plane,  Eq.  (3)  yields 
(for  the  quaternary  bowing)  1/16  ^43  + ^^3) 

and  Eq.  (5)  gives  1/8  (C^^  + ^^^3  + t C ),  which  is 

twice  as  large.  Thus,  at  the  center,  Eq.  (3)  gives  the 
average  of  the  bowing  contributions  from  the  four  ternaries, 
whereas  Eq.  (5)  gives  the  average  of  the  bowing  contribu- 
tions from  the  two  sublattices.  This  is  the  principle 
difference  between  the  interpolation  schemes  of  Refs.  6 
and  7 . 
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Comparison  with  Experiment 

Onton  and  Chicotka  compared  their  interpolation 
approach  (solution  to  Laplace's  equation)  with  measured 
values  of  the  quaternary  bandgap  in  Ga^  x^^x^l 
found  the  interpolated  values  to  be  within  ±30  meV  of  the 
measured  values.  We  have  repeated  this  comparison,  using 
their  measured  data.  The  results  are  given  in  Table  1. 
The  average  error  and  standard  deviation  of  the  error 
are  also  shown,  and  it  is  seen  that  Eqs . (3)  and  (5)  give 
comparable  errors. 


Table  I.  Comparison  of  Eqs  (3)  and  (5)  with  Experimental 

Data  for  Ga,  In  P,  As  . 

1-x  X 1-y  y 


Composition  Bandgap  (eV) 


X 

y 

Eq  (3) 

Eq  (5) 

Measured  [15] 

0.004 

0.260 

2.357 

2.353 

2.222 

0.006 

0.340 

2.244 

2.239 

2.124 

0.006 

0.420 

2.136 

2.131 

1.993 

0.009 

0.500 

2.027 

2.020 

1.881 

0.008 

0.590 

1.913 

1.907 

1.738 

0.010 

0.710 

1.763 

1.755 

1.592 

0.012 

0.830 

1.619 

1.609 

1.475 

0.017 

0.910 

1.521 

1.510 

1.351 

0.015 

0.970 

1.455 

1.447 

1.326 

0.060 

0.220 

2.336 

2.297 

2.203 

0.060 

0.280 

2.255 

2.214 

2.135 

0.070 

0.340 

2.160 

2.113 

2.038 

0.080 

0.440 

2.014 

1.963 

1. 869 

0.110 

0.650 

1.709 

1.650 

\ 1.535 

0.110 

0.700 

1.649 

1.591 

1.472 

0.170 

0.850 

1.383 

1.326 

1.236 

0.190 

0.97“^ 

1.184 

1.166 

1.165 

0.250 

0.140 

2.124 

2.054 

2.101 

0.230 

0.190 

2.100 

2.021 

2.041 

0.250 

0.240 

2.008 

1.921 

1.953 

0.250 

0.310 

1.924 

1.830 

1.809 

0.260 

0.630 

1.523 

1.430 

1.386 

0.  330 

0.800 

1.225 

1.145 

1.170 

0. 340 

0.090 

2.029 

1.971 

2.077 
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Table  I continued. 


Composition  Bandgap  (eV) 


X 

y 

Eq  (3) 

Eq  (5) 

Measured  [15, 

0.370 

0.120 

1.954 

1.883 

1.978 

0.410 

0.170 

1.843 

1.757 

1.850 

0.400 

0.030 

1.9  83 

1.957 

2.086 

0.460 

0.040 

1.887 

1.854 

2.019 

0.500 

0.040 

1.830 

1.797 

1.953 

0.570 

0.060 

1.723 

1.678 

1.854 

0.650 

0.060 

1.625 

1.580 

1.750 

0.710 

0.060 

1.558 

1.513 

1.686 

0.830 

0.100 

1.401 

1.346 

1.470 

0.  870 

0.130 

1.  333 

1.276 

1.389 

0.920 

0.190 

1.216 

1.166 

1.268 

Average 

error  ( eV ) 

-0.052 

-0.008 

rms  error  (eV) 

0.109 

0.118 

Calculation 

of  Lattice 

Constant  and 

Bandgap 

For  the  calculation  of  the  lattice  constant,  the 
three  methods  discussed  above  are  identical,  as  given  in 
Eq.  (4).  For  the  calculation  of  the  bandgap,  we  have 
elected  to  use  Moon's  procedure,  since  it  has  some 
theoretical  basis,  whereas  the  other  two  are  more  or  less 
ad  hoc. 

The  quaternary  lattice  constant  and  energy  bandgap 
contours  are  presented  in  Figure  l(a)-(i).  These  con- 
tours were  obtained  by  numerical  solutions  of  Eqs . (4) 
and  (5),  using  the  data  given  in  Table  2.  In  all  cases, 
the  lowest  quaternary  bandgap  is  plotted  in  Figure  1. 

The  shaded  regions  represent  compositions  for  which  the 
quaternary  alloy  is  an  indirect  bandgap  material. 
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Table  II.  Parameters  Used  in  the  Calculation  of 

Quaternary  Bandgaps  and  Lattice  Constants 

tendgaps  (ev)  Lattice 

O 

Compound  IXOOO)  X(IOO)  L(lll)  Constant  (A) 


AlP 

3.6 

2.42 

4.0 

5,462 

AlAs 

2.95 

2.153 

3.3 

5.6611 

AlSb 

2,25 

1.52 

1.85 

6.135 

GaP 

2.74 

2.  26 

3.0 

5.4495 

GaAs 

1.439 

1.961 

1.769 

5.64191 

GaSb 

0.69 

l.U 

0.765 

6,094 

InP 

1.34 

2.04 

1.87 

5.86875 

InAs 

0.359 

2.1 

1.6 

6.0584 

InSb 

0.175 

1.0 

0.63 

6.47877 

a) 

Binary  bandgaps  and  lattice  constants 

Alloy 

Bowing  Parameter 

r(ooo) 

X(IOO) 

L(lll) 

(Al,Ga)P 

0.0* 

0.0 

0.0 

(Al,Ga)As 

0.26 

0.02 

0.45 

(Al,Ga)Sb 

0.0 

0.0 

0.0 

(Al,In)P 

0.0 

0,0 

0.0 

(Al,In)As 

0.52* 

0.0 

0.0 

(Al,In)Sb 

0.42* 

0.0 

0,0 

(Ga,In)P 

0.758 

0,15 

0.68 

(Ga,In)As 

0.6 

0.15 

0.5 

(Ga,In)Sb 

0.43 

0.24 

0.33 

AKP.As) 

0.22* 

0.0 

0.0 

Al(P,Sb) 

1 . 2* 

0,0 

0.0 

AKAs.Sb) 

0.72* 

0.0 

0.0 

Ga(P,As) 

0.21 

0.21 

0.25 

Ga(P,Sb) 

1 . 2* 

",  0.0 

0.0 

Ga(As ,Sb ) 

0.6  5* 

0.0 

0.0 

In(P ,As ) 

0.27 

0.27 

0.26 

In(P,Sb) 

1.2 

0.0 

0.0 

In(As ,Sb) 

0.596 

0.6 

0.55 

b)  Ternary  bowing  parameters 


"'denotes  estimate  from  Figure  2.  In  the  X&L  valleys 
unknown  parameters  are  equated  to  zero. 
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Relatively  few  of  the  III-V  ternary  bowing  constants 
are  known  with  any  certainty,  especially  for  the  X(IOO) 
and  L(lll)  valleys.  Figure  2 shows  the  reported  bowing 
parameters  for  several  ternary  materials  as  a function 
of  the  lattice  constant  difference  between  the  two  end- 
point binary  compounds.  While  there  is  some  uncertainty 
in  the  experimental  data  there  does  appear  to  be  a 
definite  trend  toward  larger  bowing  parameters  with 
larger  lattice  constant  differences.  The  solid  line 
relationship  shown  in  Figure  2 has  been  used  to  estimate 
several  unknown  bowing  constants,  as  identified  by  the 
asterisks  in  Table  2.  Since  bowing  is  known  to  occur, 
it  was  felt  that  this  procedure  is  better  than  arbitrarily 
setting  the  unknown  bowing  constants  to  zero. 


I- 

[; 


Summary 


Calculated  bandgap  and  lattice  jonstant  contours  have 
been  presented  for  nine  quaternary  III-V  material  systems. 
The  calculated  values  are  based  upon  interpolation  tech- 
niques which  use  known  values  of  ternary  III-V  parameters 
to  estimate  the  bandgap  and  lattice  constant  parameters  for 
the  quaternary  systems.  Since  there  is  some  uncertainty  in 
both  the  ternary  parameters  and  the  best  interpolation 
technique  these  calculated  values  must  be  considered  as 
first  order  approximations  until  more  experimental  data 
is  obtained  for  the  quaternary  systems.  However,  tVie 
curves  should  prove  useful  for  many  studies  of  the  general 
properties  of  the  increasingly  important  III-V  quaternary 
materials . 
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BOWING  PARAMETER  (eV) 


0 8 


Figure  2. 


Experimental  bowing  parameters  versus  lattice 
constant  difference  for  III-V  ternary  alloys, 
(a)  Al,  Ga  As,  (b)  Al,  Ga  Sb , (c)  Ga,  In  P 

1-X  X 1-X  X 1-X  X 

(d)  Ga,  In  As,  (e)  Ga,  In  Sb , (f)  GaP,  As 
1-X  X 1-X  X 1-X  X 


(g)  InP,  As  , (h)  InAs,  Sb  . 
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:iiou  :-itLD  and  i;ici'.u'.;ave  lUDiLr  ii i'.ioij,Aii(.E 

IN  TERNARY  AND  r'JAl  I I'llAUY  III-V  ‘^EMICONI  UCTORS 

H.  A.  I.lttlcjohn,  J.  R.  Riuoor,  T.  H.  GUsr-on  and  L.  A.  Arlrd'jo 
f'orth  Carolina  State  University 
Raleigh,  N.  C.  27C07 


INTRODUCTION 


The  Ga7\s  KESFET  has  become  widely  recognized  as  an  excellent  micro- 
wave transistor  for  both  low-noise,  srall  signal  amplification  and  power 
amplification  through  X-band  [1,2],  arid  it  is  also  e.-.erging  as  an 
important  logic  element  fer  high-speed  digital  applications  [3].  The 
advantages  of  this  device  can  be  attributed  primarily  to  its  tcchnolcgi- 
cally  simple  structure  and  to  the  r.  :tcrial  transport  properties  of  GaAs 
[4].  Th'S  naturally  prompts  an  examination  of  other  materials  having 
properties  r.'hich  could  lead  to  improvements  in  basic  device  performance. 
Based  on  this  consideration  microwave  transistors  have  been  fabricated 
from  InP  [5]  and  Ga^_^In^As  [7]  and  theoretical  calculations  have  sug- 
gested that  InP,  As  [7]  end  Ga,  In  P,  As  [8]  have  transport  properties 
which  are  favorable  for  improvements  in  r.ESFET  pc-rformance. 

This  paper  will  discuss  KiSFET  performance  in  the  three  quaternary 

material  systems  Cc,  ^In^P,  _As  , Al,  In  P,  As  , and  Ga,  ^In  P,  Sb  . 

^ 1-x  X 1-y  y’  1-x  X 1-y  y 1-x  x 1-y  y 

This  discussion  v.'ill  center  around  results  obtained  from  the  f’onte  Carlo 
method.  These  results  have  ' cen  applied  in  device  r.odels  to  obtain 
performance  parameters  relative  to  GaAs.  Several  aspc  iS  related  to  the 
device  performance  will  also  be  addressed.  These  include  the  influence  of 
alloy  scattering,  non-unifonr  field  distributions,  non-equilibrium  tran- 
sicn.  effects  (velocity  overs: oot),  end  the  revised  conduction  band 
structure  of  GaAs  [9], 


r.MERIAL  PROPERTIES  FOR  I’.ONTE  CARLO  Of  CL'LATIONS 


If  the  I'.onte  C rlo  method  is  to  be  used  as  a sredictive  tool  for 
ternary  and  cuaternary  material:,  an  iiroortent  factor  relates  to  the 
r. liner  in  v..  ich  -aterial  param-eers  are  cstinted  and  cc-puted.  This 
can  be  illustrated  by  rcfcrcnc  to  Figure  1,  \.:i:re  the  cernosi tional 


V 


plane  (0  c . yl  , 0 <p’  <1 ) is  s!’.e.  n for  a cuet-rntry  alloy  Ai  y^y‘ 


u 

X 

First,  it  ■■  s assumed  thet  all  I ’.tcrial  p:ram:ters  n.Tced  for  the  ent- 
Carlo  r.ethtt  are  Lnewn  for  tne  binary  Ill-V  co  ipcunds  at  tlic  corn'TS 
of  the  ccmpc’Siti  -nal  plc.ne  (A-..,i:'..,AD,uD).  This  is  a reason t. ble  assuv-P-ion, 
although  some  material  p:  rar.eters  for  the  higher  energy  conduction  band 
minire  (c.c.  t a L and  X conduction  bands)  are  uncertain.  These 
parc-"t:rs  are  usually  esti  ted  free  those  for  ret'-rnls  having  these 
rimr  . a:  . lowst  C'  .cticn  ‘ nd  valley  '10,11].  > t , the  i it’- rial 

properties  rur  tne  fr  - t rn  ^y  t f al  s c I'lrul,;''"'!  t JS’  d un 


If.  ort  t 1' 

fv  " : ■ 


■Is 


: .1 3]  u-  1 
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The  final  step  In  the  procedure  Is  to  use  an  Interpolation  equation  for 
the  quaternary  material  parameter  at  the  general  composition  point  (x,y), 
or  else  to  use  available  theoretical  results  [8],  The  measure  of  success 


of  this  procedure  must  depend  on  the  accumulation  of  experimental  data, 
and  the  results  for  the  Ga,  „In„P,  As  quaternary  system  support  the 


and  the  results  for  the  Ga^_^In^P^ _^s^  quaternary  system  support  the 
validity  of  this  Interpolation  method.  These  results  Include  the  energy 


band  gap  data  [14],  the  LO  phonon  data  [15],  and  the  effective  mass. 


low-field  mobility,  and  drift-velocity  data  [16], 


DEVICE  MODELS  AfJD  MONTE  CARLO  CALCULATIONS 

V 


The  availability  of  material  parameters  allows  a steady-state  drift 
velocity-electric  field  characteristic  to  be  computed  for  a given 
material  by  the  Monte  Carlo  method.  The  relationship  of  a computed 
velocity-field  characteristic  to  two  other  characteristics  which  have 
been  used  in  device  models  to  calculate  MESFET  performance  parameters  is 
shown  In  Figure  2 for  GaAs.  The  Monte  Carlo  velocity-field  curve  has 
been  computed  using  the  r-L-X  conduction  band  ordering  for  GaAs  [9], 
and  the  material  parameters  have  been  adjusted  to  give  a good  fit  to 
experimental  velocity-field  data  for  GaAs  [21].  The  peak  velocity  is 
somewhat  lower  than  previous  calculations,  and  has  some  slight  effects 
on  the  results  presented  here.  The  device  models  have  been  described 
by  Lehovec  and  Zuleeg  [17]  and  by  Pucel  , £t  [1,18],  and  both  have 

been  used  recently  in  the  literature  as  design  models  for  MESFET's 
[1,18,19].  While  these  models  include  velocity  saturation  in  two 
different  ways,  they  do  not  include  the  negative  differential  mobility 
present  in  the  static  velocity-field  curve.  These  models  are  used  here 
in  only  a quantitative  way  to  obtain  a comparison  between  the  perforrr.ance 
parameters  for  GaAs  and  the  quaternary  materials. 

Table  1 shows  a comparison  between  the  calculated  performance 
parameters  for  a GaAs  MESFET  using  the  Monte  Carlo  velocity-field  curve 
to  provide  the  data  for  the  model  characteristics.  The  most  serious 
discrepancy  between  the  two  models  lies  in  the  values  of  the  small 
signal  drain-source  resistance,  r^^.  This  discrepancy  affects  other 

parameters  leading  primarily  to  differences  between  f and  unilateral 

gain  for  the  two  models.  This  points  out  one  of  the  most  serious 
limitations  of  these  analytical  device  models. 


TERNARY  AND  QUATERNARY  MESFET  PERFORMANCE 


The  velocity-field  characteristics  for  the  quaternary  materials 
^l-x^"x^l -y^V’  ^^l-x^"x^l-y^V  ^^1 -x^"x^l -y^^y  studied 

in  this  work.  For  the  purpose  of  this  paper  Table  2 summarizes  the  re- 
sults for  six  significant  compositions  from  these  materials.  Tabulated 
for  comparison  are  the  results  for  GaAs.  The  table  shows  velocity-field 
data  both  with  and  without  alloy  scattering  to  illustrate  the  effects  of 
this  important  scattering  process  on  the  transport  characteristics  [20]. 
y re  will  t)e  said  on  this  poirt  in  *he  ru  . i s.  ^vion. 
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Table  3 shows  the  results  of  the  MESFET  calculations  using  the 
material  parameters  In  Table  2 as  a basis.  Results  are  given  for  Pucel's 
model  only,  both  with  and  without  alloy  scattering.  There  are  slight 
discrepancies  between  the  two  models  used  for  Table  1,  but  the  results 
shown  here  are  representative  of  both  models.  It  can  be  seen  that  there 
are  predicted  Improvements  In  MESFET  performance  for  these  ternary  and 
quaternary  materials  with  respect  to  GaAs. 


OTHER  DEVICE-RELATED  ASPECTS 


There  are  several  other  aspects  of  the  Monte  Carlo  method  which  are 
viewed  as  Important  In  their  relation  to  device  performance.  A brief 
discussion  of  these  follows: 

a.  Alloy  Scattering:  The  results  have  Indicated  that  alloy  scattering 
degrades  device  performance.  Presently,  there  is  a great  deal  of  uncer- 
tainty In  the  scattering  potential  used  to  estimate  the  influence  of  alloy 
scattering  on  transport  properties,  and  three  different  scattering  pot^  - 
tials  have  been  proposed  [20].  Figure  3 shows  the  effect  of  these  scat- 
tering potentials  on  the  velocity-field  calculations  for  A1  25!^  75P  25  ~ 

As  y^and  Table  4 shows  the  extent  of  device  performance  degradation  for 

this  quaternary  material.  Depending  on  the  scattering  potential  used, 
the  degree  of  degradation  can  be  quite  severe.  The  actual  extent  of  alloy 
scattering  in  the  ternary  and  quaternary  materials  will  most  likely  deter- 
mine their  usefulness  for  theMESFET  and  other  microwave  devices. 

Transient  Transport  Effects:  The  dynamics  of  electron  transport  in 
small  devices  such  as  the  MESFET  have  been  proposed  as  a basic  frequency- 
determining phenomena  in  these  devices  [22,23].  These  effects  have  been 
investigated  in  the  quaternary  systems,  and  a typical  result  is  shown  in 

Figure  4.  Here  the  "velocity  overshoot"  is  shown  for  Ga  ^As  g 

with  and  without  alloy  scattering  and  compared  to  GaAs.  The  results  for* 
GaAs  are  slightly  different  from  previous  results  [22,23]..  The  differences 
can  be  attributed  .to  the  use  here  of  a three-band  model  appropriate  to  the 
band  structure  of  Aspnes  [9].  Again,  the  influence  of  alloy  scattering  is 
evident. 

c.  Non-Uniform  Field  Distribution:  In  a realistic  device,  the  electric 

field  intensity  will  vary  with  position  along  the  channel.  We  have 

developed  a Monte  Carlo  program  which  allows  calculations  to  be  made  with 
a non-uniform  field.  Figure  5 shows  the  velocity  along  the  channel  of  a 
device  with  a field  distribution  which  varies  linearly  from  1 kV/am  at 
x=0  to  16  kV/cm  at  x=1.5  pm  and  x=3  pm.  Also,  shown  in  this  figure  is  the 
static  drift  velocity  along  the  channel.  It  can  be  seen  that  there  are 
velocity  overshoot  effects  for  this  non-uniform  field  variation  where 
the  electron  has  been  injected  at  x=0  from  a static  electron  distribution 
corresponding  to  a field  of  1 kV/cm.  These  overshoot  effects  are  not 
as  dramatic  as  those  for  a step-field  distribution  [22,23],  but  when 
averaged  over  the  channel  could  lead  to  reduced  transit  time. 

The  overshoot  effects  for  a device  with  an  electric  field  distribution 
calculated  from  a short  channel  JFET  model  [24]  are  shown  In  Figure  6 • 

This  field  distribution  allows  for  both  source  and  drain  parasitic  regions 
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to  be  included.  For  this  model,  the  gain-bandwidth  product  has  been  numer- 
ically evaluated  and  the  results  are  indicated  in  Figure  5 for  both  GaAs 
and  Ga  73P  4^5  g.  The  results  agree  reasonably  well  with  the  device 

models,  and  are  somewhat  less  than  previous  overshoot  calculations  [25]. 


SUMMARY 


; . High  field  transport  properties  of  several  ternary  and  quaternary 
III-V  semiconductors  used  in  small  signal  device  models  predict  that 
improvements  in  MESFET  performance  over  GaAs  can  be  achieved.  The  effect 
of  alloy  scattering  generally  degrades  device  performance  and  limits 
the  improvements  which  are  predicted  for  these  materials.  In  addition, 
dynamic  effects  in  these  ternary  and  quaternary  III-V  semiconductors  are 
enhanced  when  compared  to  similar  effects  for  GaAs. 
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TABLE  1;  Device  performance  parameters  for  GaAs  obtained  by  using  Monte 
Carlo  data  in  the  models  of  Lehovec  and  Zuleeg  [17]  and  Pucel, 
etal[18]. 


A.  Channel 

width  = 

0.3wm,  Channel 

length  = 

1 .5pm, 

Channel  depth 

= 300;in, 

Np  * 10  ^ 

> 

Cg5(pf)  rjj(ohms) 

ft(GHz) 

U(db) 

la.  30.5 

.364 

445 

13.3 

24.6 

9.25 

14.9 

lb.  25.6 

.37 

460 

11.0 

18.9 

11.2 

12.7 

2a.  37.8 

.369 

1583 

16.3 

63.1 

6.7 

23.7 

2b.  31.8 

.38 

2538 

13.2 

59.2 

8.5 

23.0 

B.  Channel 

width  ■= 

0.2pm,  Channel 

length  = 

O.ZSpm, 

, Channel  depth  “ 300pra 

Nd  = 10 

cm  . 

la.  33.6 

.192 

460 

27.9 

54.8 

4.8 

21.6 

lb.  28.1 

.194 

474 

23.0 

42.0 

5.7 

19.3 

2a.  45.0 

.240 

872 

29.9 

93.6 

3.6 

27.2 

2b.  37.0 

.245 

1422 

24.1 

87.5 

4.6 

26.4 

Notation: 

g = transconductance, 
■^m 

'gs~  93te-source  capacitance. 

'•ds= 

signal  drain-source  resistance,  f^=  gain-bandwidth  product, 


f ,^=  maximum  frequency  of  oscillation,  T*=  source-drain  transit 

time,  U = Unilateral  gain  at  5GHz . Also,  gate  voltage  = 0 volts 

and  drain-source  voltage  equals  the  pinch-off  voltage.  Here 

the  number  1 denotes  Lehovec  and  Zuleeg  and  number  2 denotes 

Pucel,  et  al . The  letter  a denotes  v , = v . and  the  letter 
sat  peak 

b denotes  v^^^  = 0.8  where  is  the  maximum  Monte  Carlo 

velocity  and  is  the  device  model  saturation  velocity. 
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II 


TABLE  4:  Device  performance  parameters  for  A1  25^”  25^^  75 

illustrating  the  effects  of  different  alloy  scattering 
potentials.  The  calculations  use  the  model  of  Lehovec  and 
Zuleeg  [17]  and  are  for  a device  with  a 1.5  pm  gate  length 

17  -3 

and  channel  doping  =10  cm  . The  GaAs  parameters  are  shown 
for  comparison. 


g„(">s) 

Cg3(pf) 

r^^(ohms) 

ft (GHz)  f 

T^(psec 

a)  No  alloy 

57.9 

.407 

235 

22.6 

41.7 

5.34 

b)  AUg^ 

56.7 

.408 

236 

22.1 

40.5 

5.46 

c)  AU„ 

47.5 

.406 

294 

18.6 

34.8 

6.49 

d) 

42.7 

.410 

303 

16.6 

29.8 

7.28 

e)  GaAs 

30.5 

.364 

445 

13.3 

24.6 

9.25 

Notation: 

(see  ref. 

20):  AU^^ 

- scattering 

potential 

using  electron- 

egativity  difference,  - using  electron  affinity  difference. 


aU^q-  using  band  gap  difference. 
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Figure  1.  Quaternary  composition- 
al plane  illustrating  inter- 
polation procedure.  The  Tij'a 
are  ternary  parameters  and  Q 
is  the  interpolated  qxiater- 
nary  parameter. 
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Figure  2.  Velocity-field 
characteristic  from 
Monte  Carlo  method 
(a)  compared  to  app-X)- 
ximations  for  device 
models  (B  4 C) . 
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[ 

I Energy  band  gap  and  lattice  constant  contours  are  presented 

for  the  nine  quaternary  alloys  formed  from  Al,  Ga,  In  and 
P,  As,  Sb.  The  quaternary  bandgaps  were  obtained  using  an 
interpolation  formula  proposed  by  Moon  a^.  The  quater- 
nary lattice  constants  were  obtained  by  use  of  a linear 
interpolation  technique  using  the  binary  lattice  constants 
as  boundary  values. 

Key  words:  quaternary  alloys,  bandgap,  lattice  constant. 

Introduction 

There  has  been  considerable  interest  in  the  quater-  ^ 

nary  III-V  semiconductor  materials  for  many  applications, 

, , such  as  electro-optics  and  microwave  devices  [1-6].  One 

\ ’ significant  reason  for  this  intei-'est  is  the  ability  to 

5 synthesize  materials  with  a fixed  lattice  constant  and  a 

5 variable  range  of  energy  bandgaps  or,  conversely,  to 

? synthesize  materials  with  a constant  energy  bandgap  and 

a variable  lattice  constant  [7]. 

I *This  work  was  supported  by  a research  grant  from  the 

: Office  of  Naval  Research,  Washington,  D.C. 
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In  order  to  predict  the  energy  bandgap  or  lattice 
constant  for  a quaternary  material  a tvio-step  sequence  is 
involved.  Initially,  it  is  necessary  to  have  available 
experimental  (or  theoretical)  values  of  these  material 
parameters  for  the  four  possible  binary  III-V  constituents 
of  the  quaternary  compound.  These  values  are  readily 
available  for  most  III-V  materials  in  the  literature  [8]. 
There  are  also  four  possible  terncuny  (pseudobinary)  III-V 
systems  in  a given  quaternary  compound.  In  the  first  step, 
the  energy  bandgap  and  lattice  constant  as  a function  of 
alloy  composition  for  each  ternary  system  is  computed 
from  the  known  binary  material  parameters.  This  computa- 
tion is  on  a good  experimental  and  theoretical  basis 
[9-13]  involving  the  use  of  Vegard’s  law  and  the  concept 
of  the  "bowing  parameter"  for  the  energy  bandgap  [14]. 

The  final  step  involves  the  use  of  the  compositional 
dependence  of  the  energy  bandgap  and  lattice  constant  for 
the  four  ternary  combinations  to  interpolate  to  the  desired 
quaternary  energy  bandgap  and  lattice  constant.  This  final 
interpolation  step  is  somewhat  empirical  in  nature.  There 
appear  to  be  several  such  interpolation  techniques  used  in 
the  literature  [6,7,15]  and  while  these  are  similar,  they 
also  have  non-trivial  differences.  This  paper  will  discuss 
this  critical  interpolation  procedure  in  view  of  the  small 
amount  of  experimental  quaternary  data,  and  will  present 
calculations  for  energy  bandgaps  and  lattice  constants 
for  the  nine  quaternary  alloys  formed  from  Al,  Ga,  In  and 
P,  As,  Sb  using  one  of  these  interpolation  techniques. 


Notation 


There  appears  to  be  no  generally-accepted  scheme  for 

symbolically  describing  the  III-V  quaternary  alloys.  In 

tliis.  paper  the  following  convention  has  been  adopted. 

Within  both  the  group  III  and  group;  V pairing,  the  first 

position  is  occupied  by  the  element  with  lowest  atomic 

numi  er.  The  compios i t ion  variables  are  associated  with  the 

four  elements  in  the  order  1-x,  x,  1-y,  y.  For  example,  a 

.quaternary  alloy  will  be  denoted  as  A,  B C,  D . Here  A 

1-xxl-yy 

and  B are  group  III  elements  with  A having  lower  atomic 
number  than  B.  Likewise,  C and  D are  group  V elements  with 
C having  lower  atomic  numlier  than  r< . Following  previous 
auttiof.  (6,V  1,  a quaternary  a 1 loy  parameter  (e.g.,  bandgap 
or  lattice  coir.tant  ) is  descrit'ed  by  a •■urface  0(x,y)  over 
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In  order  to  predict  the  energy  bandgap  or  lattice 
constant  for  a quaternary  material  a two-step  sequence  is 
involved.  Initially,  it  is  necessary  to  have  available 
experimental  (or  theoretical)  values  cf  these  material 
parameters  for  the  four  possible  binary  III-V  constituents 
of  the  quaternary  compound.  These  values  are  readily 
available  for  most  III-V  materials  in  the  literature  [8]. 
There  are  also  four  possible  ternary  (pseudobinary)  III-V 
systems  in  a given  quaternary  compound.  In  the  first  step, 
the  energy  bandgap  and  lattice  constant  as  a function  of 
alloy  composition  for  each  ternary  system  is  computed 
from  the  known  binary  material  parameters.  This  computa- 
tion is  on  a good  experimental  and  theoretical  basis 
[9-13]  involving  the  use  of  Vegard's  law  and  the  concept 
of  the  "bowing  parameter"  for  the  energy  bandgap  [14]. 

The  final  step  involves  the  use  of  the  compositional 
dependence  of  the  energy  bandgap  and  lattice  constant  for 
the  four  ternary  combinations  to  interpolate  to  the  desired 
quaternary  energy  bandgap  and  lattice  constant.  This  final 
interpolation  step  is  somewhat  empirical  in  nature.  There 
appear  to  be  several  such  interpolation  techniques  used  in 
the  literature  [6,7,15]  and  while  these  are  similar,  they 
also  have  non-trivial  differences.  This  paper  will  discuss 
this  critical  interpolation  procedure  in  view  of  the  small 
amount  of  experimental  quaternary  data,  and  will  present 
i calculations  for  energy  bandgaps  and  lattice  constants 

for  the  nine  quaternary  alloys  formed  from  Al,  Ga , In  and 
P,  As,  Sb  using  one  of  these  interpolation  techniques. 

Notation 

There  appears  to  be  no  generally-accepted  scheme  for 
symbolically  describing  the  III-V  quaternary  alloys.  In 
this. paper  the  following  convention  has  been  adopted. 

Within  both  the  group  III  and  group  V pairing,  the  first 
position  is  occupied  by  the  element  with  lowest  atomic 
numi:'er.  The  composition  variables  are  associated  with  the 
four  elements  in  the  order  1-x,  x,  1-y , y.  For  example,  a 
.quaternary  alloy  will  be  denoted  as  x^x^l  ^ 

and  B are  group  III  elements  with  A having  lower  atomic 
number  than  B.  Likewise,  C and  D are  group  V elements  with 
C liaving  lower  atomic  num!ier  than  li . Following  previous 
authors  [6,7],  a quaternary  alloy  paramrter  (e.g.,  bandgap 
or  lattice  constant)  is  described  by  a '.urface  0(x,y)  over 
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the  x,y  composition  plane  (0<x<l,  0<y<l).  At  the 
corners  (x  and  y equal  to  zero  and/or  one)  the  values  of 
the  parameter  for  the  four  binary  elements  Q(0,0)  = B. , 
Q(l,0)  = B , Q(l,l)  = and  Q(0,1)  = are  obtained. 
Along  the  Boundaries  of  the  plane  the  parameter  for  the 
four  ternary  elements  Q(x,0)  = T (x),  Q(l,y)  = T^,(y), 

43^ 


Q(x,l)  = Tj^q(x)  and  Q(0,yj  = T ^^y)  are  obtained. 


Estimation  of  Alloy  Parameters 

Many  ternary  alloy  parameters  (e.g.  lattice  constant) 
are  obtained  by  linear  interpolation  from  those  of  the 
constituent  binary  compounds,  i.e., 

T.  .(x)  = xB.  + (l-x)B.  (1) 

1]  1 1 

For  some  parameters  the  theoretical  variation  with  compo- 
sition is  nonlinear  [14,16].  For  the  direct  bandgap  for 
example,  Thompson  and  Woolley  [9]  have  shown  that 


T..(x)  = xB.  + (l-x)B.  - C.  .x(l-x)  (2) 

11  1 1 ID 

where  . is  the  bowing  parameter  for  the  ternary  alloy 
bandgap 

In  the  absence  of  definitive  theories  for  quaternary 
parameters,  estimates  of  a quaternary  alloy  parameter 
Q(x,y)  must  be  obtained  by  interpolation  from  the  four 
ternary  alloy  parameters  T.  . Various  interpolation 
schemes  have  been  proposed^^  Onton  and  Chicotka  [15] 
used  the  solution  of  Laplace's  equation  subject  to  the 
boundary  conditions  Q(x,0)  = T^.,(x),  etc.  In  Monte  Carlo 
transport  studies  [6]  the  present  authors  have  used  the 
interpolation  equation 


Q(x,y)  = 


x(l-x)[  ( l-v)Tj^j(x  l+yTj^gCx)  ]+y ( 1-y  )[  ( 1-x )T^i^ (y  ) + x y ) ] 
x( l-x)+y( 1-y ) 


(3) 


This  interpolation  equation  reduces  to  the  ternary  param- 
eters on  the  quaternary  plane  boundaries  and  to  the  average 
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of  the  ternary  parameters  at  the  midpoint  (x=0. 5 ,y  = 0. 5) 
of  the  compositional  plane.  This  is  incidentally  the 
solution  to  Laplace's  equation  if  the  T. . are  linear  as 
in  Eq.  (1).  Moon  [7]  have  propo^^d  similar  schemes 

for  the  estimation  of  lattice  constant  and  bandgap.  For 
the  lattice  constant,  Moon  gives  (in  the  above  notation) 

Q(x,y)  = B^+(B2-Bj^)x  + (B^-B^)y  + (B^-B2tB2-B^)xy  (4) 


For  the  lattice  constant,  the  T. . are  assumed  to  be 
linear  in  their  arguments,  so  the  solution  to  Laplace's 
equation  is  identical  to  Eq‘.  (3),  which  in  turn  also 
reduces  to  Eq.  (4). 

For  the  bandgap.  Moon  e_^  a^.  have  used  the  equation 
Q(x,y)  = (l-x)T^j^(y)  + x ^ 


where  the  T.  . are  determined  from  Eq.  (1)  and 
ID 


A = x(l-x)[(l-y)C^2'''y^43^  y ( 1-y )[ ( (6) 


An  important  difference  between  this  bandgap  estimate  and 
that  obtained  from  Eq,  (3)  is  in  the  manner  in  which  the 
ternary  bowing  parameters  C..  enter  the  calculation.  In 
Eq.  (3),  the  ternary  bowing^Jf fects  are  included  in  the 
T. .,  whereas  in  Moon's  method  the  bowing  enters  as  a 
separate  quaternary  bowing  parameter  term  given  in  Eq.  (6). 
At  the  center  of  the  composition  plane,  Eq.  (3)  yields 

: "14  * "23> 

) , which  is 


(for  the  quaternary  bowing)  1/16  (C^^  ^43 


and  Eq.  (5)  give: 


i/e  c,3 

Thu 


+ C 


14 


+ C, 


23 

twice  as  large.  Thus,  at  the  center,  Eq . (3)  gives  the 
average  of  the  bowing  contributions  from  the  four  ternaries, 
whereas  Eq . (5)  gives  the  average  of  the  bowing  contribu- 
tions from  the  two  sublattices.  This  is  the  principle 
difference  between  the  interpolation  schemes  of  Refs.  6 
and  7 . 
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Comparison  with  Experiment 


Onton  and  Chicotka  compared  their  interpolation 

approach  (solution  to  Laplace's  equation)  with  measured 

values  of  the  quaternary  bandgap  in  Ga,  In  P,  As  and 

1-x  X 1-y  V 

found  the  interpolated  values  to  be  within  ±30  mev  of  the 
measured  values.  We  have  repeated  this  comparison,  using 
their  measured  data.  The  results  are  given  in  Table  1. 
The  average  error  and  standard  deviation  of  the  error 
are  also  shown,  and  it  is  seen  that  Eqs . (3)  and  (5)  give 
comparable  errors. 


Table  I. 

Comparison 

Data 

of  Eqs  (3) 

for  Ga, 

1-x 

and  (5)  w 

In  P,  As 

X 1-y  y 

ith  Experimental 

Composition 

Ba-dgap 

(eV) 

X 

y 

Eq  (3) 

Eq  (5) 

Measured  [15] 

0.004 

0.260 

2.357 

2.353 

2.222 

0.006 

0.340 

2.244 

2.239 

2.124 

0.006 

0.420 

2.136 

2.131 

1.993 

0.009 

0.500 

2.027 

2.020 

1.881 

0.008 

0.590 

1.913 

1.907 

1.738 

0.010 

0.710 

1.763 

1.755 

1.592 

0.012 

0.830 

1.619 

1.609 

1.475 

0.017 

0.910 

1.521 

1.510 

1.351 

0.015 

0.970 

1.455 

1.447 

1.326 

0.060 

0.220 

2.336 

2.297 

2.203 

0.060 

0.280 

2.255 

2.714 

2.135 

0.070 

0.340 

2.160 

2.113 

2.038 

0.080 

0.440 

2.014 

1.963 

1.869 

0.110 

0.650 

1.709 

1.650 

\ 1.535 

0.110 

0.700 

1.649 

1.591 

1.472 

0.170 

0.850 

1.383 

1.3  26 

1.236 

0.190 

0.975 

1.184 

1.166 

1.16  5 

0.250 

0.140 

2.124 

2.054 

2.101 

0.230 

0.190 

2.100 

2.021 

2.041 

0.250 

0.240 

2.008 

1.921 

1.953 

0.250 

0.310 

1.924 

1.830 

1.809 

0.260 

0.630 

1.523 

1.430 

1.388 

0.330 

0.800 

1.225 

1.145 

1.170 

0. 340 

0.090 

2.029 

1.971 

2.077 
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Table  I continued. 


Composition  Bandgap  (eV) 


X 

y 

Eq  (3) 

Eq  (5) 

Measured  [15] 

0.370 

0.120 

1.954 

1.883 

1.978 

0.410 

0.170 

1.843 

1.757 

1.850 

0.400 

0.030 

1.983 

1.957 

2.086 

0.460 

0.040 

1.887 

1.854 

2.019 

0.500 

0.040 

1.830 

1.797 

1.953 

0.570 

0.060 

1.723 

1.678 

1,854 

0.650 

0.060 

1.625 

1.580 

1.750 

0.710 

0.060 

1.558 

1.513 

1,686 

0.830 

0.100 

1.401 

1.346 

1,470 

0.870 

0.130 

1.333 

1.276 

1,389 

0.920 

0.190 

1.216 

1.166 

1.268 

Average 

error  ( eV ) 

-0.052 

-0.008 

rms  error  (eV) 

0.109 

0.118 

Calculation  of  Lattice  Constant  and 

Bandgap 

For  the  calculation  of  the 

lattice  constant,  the 

three  methods  discussed  above  are  identical 

, as  given  in 

Eq.  (4), 

, For  the  calculation  of  the  bandgap,  we  have 

elected 

to  use  Moon's  procedure 

, since  it  has  some 

theoretical  basis,  whereas  the  other  two  are  more  or  less 

ad  hoc. 

The  quaternary  lattice  constant  and  energy  bandgap 
' , contours  are  presented  in  Figure  l(a)-(i).  These  con- 

I tours  were  obtained  by  numerical  solutions  of  Eqs . (4) 

' and  (5),  using  the  data  given  in  Table  2.  In  all  cases, 

the  lowest  quaternary  bandgap  is  plotted  in  Figure  1. 

The  shaded  regions  represent  compositions  for  which  the 
quaternary  alloy  is  an  indirect  bandgap  material. 
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Table  II.  Parameters  Used  in  the  Calculation  of 

Quaternary  Bandgaps  and  Lattice  Constants 


Bandgaps 

(ev) 

Lattice 

Compound 

IXOOO) 

X(IOO) 

L(lll) 

Constant  (A) 

AlP 

3.6 

2.42 

4.0 

5.462 

AlAs 

2.95 

2.153 

3.3 

5.6611 

AlSb 

2.25 

1.52 

1.85 

6.135 

GaP 

2.74 

2,26 

3,0 

5.4495 

GbAs 

1.439 

1.961 

1.769 

5.64191 

GaSb 

0.69 

l.U 

0.765 

6.094 

InP 

1.34 

2.04 

1.87 

5.86875 

InAs 

0.359 

2.1 

1.6 

6.0584 

InSb 

0.175 

1.0 

0.63 

6.47877 

a) 

Binary  bandgaps  and  lattice  constants 

Alloy 

Bowing  Parameter 

r(ooo) 

X(IOO) 

L(lll) 

(Al,Ga)P 

0.0* 

0.0 

0.0 

(Al,Ga)As 

0.26 

0.02 

0.45 

(Al,Ga)Sb 

0.0 

0.0 

0,0 

(Al,In)P 

0.0 

0.0 

0.0 

(Al,In)As 

0.52* 

0.0 

0.0 

(Al,In)Sb 

0.42* 

0.0 

0.0 

(Ga,In)P 

0.758 

0.15 

0,68 

(Ga,In)As 

0.6 

0.15 

0.5 

(Ga,In)Sb 

0.43 

0.24 

0.33 

AKP.As) 

0.22* 

0.0 

0.0 

AKP.Sb) 

1 . 2* 

0.0 

0.0 

AKAs.Sb) 

0.72* 

0.0 

0.0 

Ga(P,As) 

0.21 

0.21 

0.25 

Ga(P,Sb) 

1.2* 

% 0.0 

0.0 

Ga(As,Sb) 

0.65* 

0.0 

0.0 

ln(P,As) 

0.27 

0.27 

0.26 

In(F,Sb) 

1.2 

0.0 

0.0 

In{ As  ,Sb) 

0.596 

0.6 

0.55 

b)  Ternary  bowing  parameters 


'■‘denotes  estimate  from  Figure  2.  In  the  X&L  valleys 
unknown  parameters  are  equated  to  zero. 
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Estimation  of  Unknown  Ternary  Bowing  Constants 

Relatively  few  of  the  III-V  ternary  bowing  constants 
are  known  with  any  certainty,  especially  for  the  X(IOO) 
and  L(lll)  valleys.  Figure  2 shows  the  reported  bowing 
parameters  for  several  ternary  materials  as  a function 
of  the  lattice  constant  difference  between  the  two  end- 
point binary  compounds.  While  there  is  some  uncertainty 
in  the  experimental  data  there  does  appear  to  be  a 
definite  trend  toward  larger  bowing  parameters  with 
larger  lattice  constant  differences.  The  solid  line 
relationship  shown  in  Figure  2 has  been  used  to  estimate 
several  unknown  bowing  constants,  as  identified  by  the 
asterisks  in  Table  2.  Since  bowing  is  known  to  occur, 
it  was  felt  that  this  procedure  is  better  than  arbitrarily 
setting  the  unknown  bowing  constants  to  zero. 


Summarv 


Calculated  bandgap  and  lattice  constant  contours  have 
been  presented  for  nine  quaternary  III-V  material  systems. 
The  calculated  values  are  based  upon  interpolation  tech- 
niques which  use  known  values  of  ternary  III-V  parameters 
to  estimate  the  bandgap  and  lattice  constant  parameters  for 
the  quaternary  systems.  Since  there  is  some  uncertainty  in 
both  the  ternary  parameters  and  the  best  interpolation 
technique  these  calculated  values  must  be  considered  as 
first  order  approximations  until  more  experimental  data 
is  obtained  for  the  quaternary  systems.  However,  the 
curves  should  prove  useful  for  many  studies  of  the  general 
properties  of  the  increasingly  important  III-V  quaternary 
materials . 
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Figure  2.  Experimental  bowing  parameters  versus  lattice 
constant  difference  for  III-V  ternary  alloys, 
(a)  Al,  Ga  As,  (b)  Al,  Ga  Sb,  (c)  Ga,  In  P, 

1-X  X 1-X  X 1-X  X 

(d)  Ga,  In  As,  (e)  Ga,  In  Sb,  (f)  GaP,  As  , 
1-X  X 1“X  X l-X  X 

(g)  InP,  As  , (h)  InAs,  Sb  . 

® i“X  X*  1— X X 
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